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ABSTRACT 



We report 4 new detections of 21-cm absorption from a systematic search of 21-cm absoiption in a sample of 17 strong (rest equivalent 
width, W r (Mg iLi2796)>l A) intervening Mg n absorbers at 0.5<z a b s <1.5. We also present 20-cm milliarcsecond scale maps of 40 
quasars having 42 intervening strong Mg n absorbers for which we have searched for 21-cm absorption. These maps are used to 
understand the dependence of 21-cm detection rate on the radio morphology of the background quasar and address the issues related 
to the covering factor of absorbing gas. Combining 21-cm absoiption measurements for 50 strong Mg n systems from our surveys with 
the measurements from literature, we obtain a sample of 85 strong Mg n absorbers at 0.5<z abs <l and l.l<z abs <1.5. We present detailed 
analysis of this 21-cm absoiption sample, taking into account the effect of the varying 21-cm optical depth sensitivity and covering 
factor associated with the different quasar sight lines. We find that the 21-cm detection rate is higher towards the quasars with flat or 
inverted spectral index at cm wavelengths. About 70% of 21-cm detections are towards the quasars with linear size, LS<100pc. The 
21-cm absoiption lines having velocity widths, AV>100km s -1 are mainly seen towards the quasars with extended radio morphology 
at arcsecond scales. However, we do not find any correlation between the integrated 21-cm optical depth, Jrdv, or the width of 21-cm 
absorption line, AV, with the LS measured from the milliarcsecond scale images. All this can be understood if the absorbing gas is 
patchy with a typical correlation length of ~30-100pc. We confirm our previous finding that the 21-cm detection rate for a given 
optical depth threshold can be increased by up to a factor 2 by imposing the following additional constraints: Mg n doublet ratio <1.1, 
W(Mg n)/W(Fe n)< 1 .47 and W(Mg i)/W(Mg n)>0.27. This suggests that the probability of detecting 21-cm absoiption is higher in 
the systems with high A'(H i). We show that within the measurement uncertainty, the 21-cm detection rate in strong Mg n systems 
is constant over 0.5<z a bs<l-5, i.e., over ~30% of the total age of universe. We show that the detection rate can be underestimated by 
up to a factor 2 if 21-cm optical depths are not corrected for the partial coverage estimated using milliarcsecond scale maps. Since 
stellar feedback processes are expected to diminish the filling factor of cold neutral medium over 0.5< z <1, this lack of evolution 
in the 21-cm detection rate in strong Mg n absorbers is intriguing. Large blind surveys of 21-cm absorption lines with the upcoming 
Square Kilometre Array pathfinders will provide a complete view of the evolution of cold gas in galaxies and shed light on the nature 
of Mg ii systems and DLAs, and their relationship with stellar feedback processes. 

Key words, quasars: absoiption lines - galaxies: evolution - galaxies: high redshift - galaxies: ISM - galaxies: star formation. 



1. Introduction 

The diffuse interstellar medium (ISM) exhibits a wide range of 
physical conditions such as temperature, density and radiation 
field, and structures at tens of AU to kpc scales in the form of 
| shells, filaments and spurs. The conditions in the ISM are inti- 
$_i , mately related to the local star-formation through various stel- 
lar feedback mechanisms. Therefore understanding the physical 
conditions and the structure of gas in the ISM, and the processes 
that maintain these is a key topic in the field of galaxy evolu- 
tion. The absorption lines seen in the spectra of distant luminous 
sources such as quasi-stellar objects (QSOs) and Gamma-ray 
bursts (GRBs) provide a powerful tool to determine the physi- 
cal and chemical state of the different phases of ISM in both the 
Galaxy and external galaxies. These absorbers provide a lumi- 
nosity unbiased view of the galaxy evolution and measurements 
of the physical parameters of gas at the scales much smaller than 
that possible with the continuum and emission line studies. 

Due to their observability with the ground based optical tele- 
scopes, the resonant absorption lines of Mg n are commonly 
used for detecting the interstellar media and/or gaseous halos 



surrounding the galaxies at z<2. The availability of large sam- 
ples of low resolution quasar spectra from the Sloan Digital 
Sky Survey (SDSS) has led to the compilation of large homo- 
geneous catalogues of Mg n absorbers at 0.35<z a b s <2.3 that are 
complete down to W r (Ms n/12796)=lA (see e.g.. iProchter et all 
120061 : lOuider et al.ll2.Ql ll) . It has been shown using the Hubble 
Space Telescope (HST) spectroscopic observations that a sig- 
nificant fraction (36%+6%) of Mg n systems with the rest 
equivalent widthsQ W r >0.5A and W r (Fe n/l2600)>0.5A have 
N(H i)>2xl0 20 cm 2 and there fore are bonafide damped Lya 
(DLA) systems dRao et alj|2006h . The optical imaging and spec- 
troscopic surveys have shown that the Mg n absorbers are asso- 
ciated with the galaxies of a wide range of morphological type 
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While it is well established that the Mg n absorption lines 
arise from the gas associated with the galaxies, the exact na- 
ture and the underlying physical process driving the gas seen 
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in absorption is still a matter of de bate. Using the integral field 
spectroscopy, iBouche et al.l (120071) have detected Ha- emission 
associated with the 14 intervening Mg n absorbers (W T > 2 A) 
within the QSO impact parameters of l-40kpc implying large 
star formation rates (SFRs) of 1-20 M yr~'. Similarly, the 
[O n] emission lines corresponding to the absorbing galaxies 
have been detected in the individual a s well as the stacked QSO 
spectra of W r >0.7 A Mg n absorbers dNoterdaeme et alj|2010bl 
iMenard et al.ll2.Ql Tb . These results point towards a strong phys- 
ical connection between the starburst phenomena and the ab- 
sorbing g as associated with the strong (W r >l A) Mg n absorbers 
(see a lso iProchter et aDl2006l iMurphv et aHl2007UNestor et al.l 
l20ll . On the other hand, the weaker Mg n ab sorbers are 
believed to be tracing the infallin g/accreting gas dChen et al.l 
1201 Obi iKacprzak & Churchillll201 ll) . 

The detection of 21 -cm absorption in a well defined sample 
of Mg ii absorbers towards the radio-loud quasars can be used 
to estimate the filling factor of cold gas in Mg n systems. Until 
blind surveys of radio absorption lines become a real possibil- 
ity, this is also the only practical way to search for the 21 -cm 
absorption at 0.2<z a b s <2 and trace the evolution of cold gas in 
galaxies. Since the cold neutral medium (CNM) filling factor of 
gas depends sensitively on the stellar feedback, such surveys can 
also provide insights into the physical processes driving the gas 
seen in Mg n absorption. There have be en several searches for 



21 -cm abso rption in Mg n systems (e.g. Briggs & Wolfdll983i 
lLandl2000i ICurran et alj|2007t iKanekar et al J2009I) . A svstem- 



atic survey of 21 -cm absorption in a sample of 35 intervening 
Mg ii abs orbers at l.l<z a hs< 1 -5 was presented in iGupta et all 
d2007l) and lGupta etal] d2009l hereafter referred to as G09). The 
detection of 9 new 21 -cm absorbers from this survey signifi- 
cantly increased the number of known 21 -cm absorbers at the 
intermediate redshifts. 

G09 provided first estimate of the number density per unit 
redshift of 21-cm absorbers («2i) associated with the W r >l A 
Mg ii absorbers. Hereafter, we will refer to W r >l A absorbers 
as the strong Mg n absorbers. Using the twenty-one 21-c m ab- 
sorptio n measurements at 0.3 <z a b s < 1 from the sample of lLand 
d2000l) . they showed that the «2i at z = 1.3 may have decreased 
with respect to the «2i at z = 0.5. The evolution was found to be 
stronger for the stronger (W r >1.8A) Mg n systems. This could 
be due to the evolution of the CNM filling factor of the gas in 
galaxies, and the fact that the stronger Mg n syst ems at higher-z 
possibly arise in the outflows d riven by starbursts (IProchter et al.l 
20061 iNoterdaeme et al.ll20l"2l) . This interpretation of the results 
from the survey was limited by the availability of the milliarc- 
second scale radio images of the background quasars and the 
lack of sensitive 21-cm absorption measurements at z < 1. 

As an extension of our work presented in G09, we present 
here the new 21-cm absorption measurements of the 17 Mg n 
systems at 0.5<z a b s <T5 and the 20-cm Very Long Baseline 
Array (VLBA) maps of the 40 quasars to study the evolution 
of CNM filling factor of Mg n systems and investigate the ef- 
fect of radio source structure on the detectability of 21-cm ab- 
sorption. The 21-cm absorption measurements were made us- 
ing the Green Bank Telescope (GBT), the Giant Metrewave 
Radio Telescope (GMRT) and the Westerbork Synthesis Radio 
Telescope (WSRT). The details of these 21-cm absorption and 
the VLBA observations, and the new 21-cm absorption detec- 
tions are provided in Section [2] In Section [3] we describe the 
assumptions and method used to determine the correction for 
partial coverage using the milliarcsecond scale VLBA maps. We 
define a homogeneous sample that is selected primarily on the 





Fig. 1. GMRT maps of J0956+4046 (left) and J2031 + 1219 that 
are resolved in our observations at ~610 MHz with contour lev- 
els as nx(-l, 1,2,4, ...) mJy beam -1 , where n=2.8 and 10 re- 
spectively. The maps have rms of 0.36 and 1.4 mJy beam -1 , 
and the beams are 5.43"x4.76" (position angle PA=37°), 
and 5.34"x4.96" (PA=-39°) respectively. The position of 
the optical source is marked as '+' (SDSS for J0956+4046; 
ICondon et alii 19771 for J2031 + 1219). 



basis of Wr >1A criterion. This sample is used to derive the 
21-cm detection rate in strong Mg n absorbers as a function 
of redshift. In Section |4] we discuss the relationship between 
the 21-cm absorption and the radio structure of the background 
quasar. The relationship of 21-cm absorption with the properties 
of metal absorption lines and the dust content of the absorbing 
gas derived from the optical spectra are discussed in Section [5] 
We discuss the implications of our results in Section [6] The re- 
sults are summarised in Section [7] 

Throughout this paper we use the ACDM cosmology with 
Q m =0.27, Q A =0.73 and H Q =71 km s -1 Mpc" 1 . 



2. Data 

2.1. Redshifted 21 -cm line observations 
2.1.1. Sample 

Our sample of Mg ii absorbers to search for 21-cm absorption 
is listed in Table Q] The sample is mostly drawn from the SDSS 
Data Release 7 using an automatic procedure that fits the QSO 
continuum redwards of the Lyo- emission using Savitzky-Golay 
filtering and identifies Mg ii doublets at 0.35<z a b s <2.3 through 
a correlation analysis. Similar to G09, we confine our sample to 
the Mg ii systems with W r >l A. The sample was designed for ob- 
servations with the GMRT 610-MHz and WSRT UHF-high re- 
ceivers (see Section l2.1.2b . Therefore we selected Mg n systems 
at 0.5<z a bs<1.0 and 1 . 1 <z a bs< 1 -5. Although Mg n systems down 
to the redshifts of ~0.35 can be detected in the SDSS spectra 
and the corresponding redshifted 2 1 -cm absorption is observable 
with the WSRT UHF-high receiver, we chose the lower redshift 
cut-off of our sample to be at z a bs=0.5 to avoid the frequencies 
around 950 MHz that are known to be affected by strong radio 
frequency interference (RFI). By cross correlating the position 
of the QSOs having strong intervening Mg n absorption with 
radio sources in the Faint Images of the Radio Sky at Twenty- 
Centimeters (FIRST) catalogue we identified Mg n systems that 
are in front of the quasars brighter than lOOmJy at 20-cm. After 
visual inspection of the optical spectra and radio images to en- 
sure that there are no false identifications, we selected 19 Mg n 
systems that are not covered in G09 for our observations. 

In addition, we selected 4 strong Mg n systems towards 
well known blazars with the metal absorption line prop- 
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Table 1. The sample of strong Mg n systems (W r (Mg n/l2796)>l A) observed with the GBT, GMRT and WSRT The systems for 
which no useful 21 -cm absorption spectra could be obtained are also given. Listed from left to right are the names of the quasars, 
their emission redshifts, redshifts of the intervening Mg n system, rest frame equivalent widths of the Mg n/12796, Mg ii/i2803, 
Mg 1/12852 and Fe n/12600 absorption lines, and 1.4 GHz peak flux densities of the quasars as determined from the FIRST or NVSS 
surveys. The last column specifies whether the q uasar is compact (C ) or resolved (R) at the FIRST resolution. The quasars with the 
deconvolved sizes <2" are classified as compact dWhite et alj|1997l) . 



Quasar 




~abs 


W r (Mg H/12796) 


W r (Mg H/12803) 


W r (Mg 1/12852) 


W r (Fe H/12600) 


Flux 


Mor. 








(A) 


(A) 


(A) 


(A) 


(mjy/beam) 










Systems with 21-cm absorption spectra 








T045701— 919459 


1 00? 


8999 


2.20+0.02 


1.90±0.02 


0.61±0.02 


1 S7+D 09 


1 640 




torooi6+50i 044 


1615 


1 4146 


1.94±0.04 


1.70±0.04 


0.45±0.05 


1 .Z.UIEU.UJ 


1119 


V 


T081 71 0+915991 


1.732 


1.3060 


1.70±0.13 


1.53±0.13 


0.80±0.15 


1.38+0.17 


210.4 


Q 


1093035+464408 


2.033 


0.6216 


3.18+0.14 


2.87±0.14 


1.21+0.14 


2.44+0.13 


313.9 


c 


J09563 1+404628 


1.510 


1.3237 


5.02±0.37 


4.19±0.35 


2.43±0.36 


1.97+0.57 


30.6 


R 


T 1007 18+225 126 


2.308 


0.5602 


1.19±0.32 


1.89±0.24 


0.64±0.27 


1.03±0.19 


341.3 


c 


11 14856+525425 


1.632 


0.8306 


2.77±0.04 


2.39±0.05 


0.76±0.06 


1.66±0.05 


95.6 


c 


1191604+584333 


1.454 


0.7487 


1.55±0.12 


1.44±0.12 


0.47±0.12 


0.97±0.11 


376.8 


R 


1125227+442737 


1.347 


0.9112 


1.39+0.15 


1.33±0.14 


0.36±0.13 


1.24±0.14 


248.7 


R 


1132901 + 105304 


1.933 


0.6715 


1.47±0.18 


1.45±0.18 


<0.24 


1.00+0.25 


111.2 


c? 


J133335+ 164903 


2.089 


0.7448 


1.80±0.04 


1.46+0.04 


0.18±0.04 


0.53±0.05 


409.1 


c 


J 140856-075227 


1.500 


1.2753 


2.08±0.05 


1.90+0.05 


0.62±0.05 


1.53±0.05 


604.0 


c 


J141030+614136 


2.247 


0.7596 


1.31±0.16 


1.48+0.16 


0.57±0.10 


0.84±0.19 


116.1 


c 


J150124+561949 


1.467 


1.2788 


1.38+0.10 


1.02±0.09 


0.37±0.10 


0.72±0.12 


148.8 


c 


J163638+211255 


1.801 


0.8000 


2.26±0.20 


2.41+0.21 


0.53±0.24 


1.62±0.20 


365.8 


c 


J203155+121941 


1.215 


1.1158 


1.29±0.05 


1.16+0.05 


0.54±0.05 


1.12+0.05 


941.0 


Rt 


J212912-153841 


3.268 


0.6628 


2.49±0.003 








589.7 


R ; 






Systems with VLBA image but no 21-cm absorption spectra 






J08 1534+330529 


2.426 


0.8515 


2.67±0.28 


2.19+0.31 


0.55±0.15 


1.12+0.27 


327.9 


C 


Jl 10021+162914 


3.380 


0.8540 


1.99+0.11 


1.64±0.11 


0.66±0.12 


1.08+0.13 


256.5 


C 


Jl 15734+ 163859 


1.061 


0.7624 


1.19+0.07 


1.05+0.07 


0.25±0.08 


0.56±0.07 


734.3 


C 


J121332+130720 


1.139 


0.7718 


1.19+0.08 


1.10+0.08 


0.36±0.09 


0.81+0.08 


1281.0 


C 


J130036+082802 


1.090 


0.8665 


1.39±0.08 


1.07±0.08 


0.48±0.10 


0.78+0.11 


101.3 


C 


J143009+104326 


1.710 


1.2431 


1.64±0.09 


1.53±0.09 


<0.10 


0.64+0.17 


303.9 


c 



Notes. (t) Not covered in FIRST. The flux densities for J0457-2324, J2031 + 1219 and J2129-1538 are taken from the NVSS. The morphol- 
ogy classification for J0457-2324 and J2129— 1538 is based on the 1.4GHz VLA images with the resolution of ~3" and ~1.4" respectively 
dUlvestad et alJI 19811 : iNeff & Hutchingslll99fJ) . More than 90% of the flux density of J0457-2324 is contained in an unresolved component 
whereas J2129— 1538 exhibits a core-jet morphology with linear size of 4.6" and 95% of the total flux density in the core component. The quasar 
J203 1 + 1219 is extended in the GMRT image at 610 MHz (resolution~5") presented in Fig.fTJ 



erties suggesting a high probability of them being DLAs. 
The z abs = 0.8922 system towards the BL Lac J0457-2324 
(PKS 0454-234), the z abs = 1.2753 system towards the FSRQ 
J1408-0752 (PKS 1406-076), and the z abs = 1.1158 system 
towards the BL Lac J2031 + 1219 (PKS 2029+121), were se- 
lec ted from the sample o f Mg n systems towards blazars given 
in iBergeron et al.l (1201 ll) . The z a bs= 0.6628 system towards 
the FSRQ J2129-1538 (PKS 212 6-158) was identified in the 
UVES/VLT large programme data (IBergeron et al.ll2004l) . 



In total, we have a sample of 23 strong Mg n systems. The 
rest equivalent widths of Mg 1/12852, Mg H/12796, Mg H/12803 
and Fe H/12600 for these are given in the Table Q] The equiva- 
lent widths have been estimated by integrating over the absorp- 
tion profile as done in G09. The Mg 1/12852, Mg ii/i2803, and 
Fe 11/12600 absorption lines for the system towards J2129-1538 
fall in the Lya forest, and their profiles are blended with the other 
absorption features. For this reason, the equivalent widths for 
these absorption lines for this system are not given in Table Q] A 
detailed analysis of this system will be presented in Boisse et al. 
(in prep.). 



2.1.2. Observations and data reduction 

The systems at 0.5<z a bs<L0 were observed using the WSRT 
UHF-high receiver whereas the systems at l.l<z a bs<L5 were 
observed using the GMRT 610-MHz receiver. The latter corre- 
sponds to the frequency coverage of the 610-MHz receiver. The 
systems at z a bs<l that could not be observed with the WSRT or 
resulted in the spectra severely affected by the RFI were sub- 
sequently observed with the GBT prime-focus PFl-800MHz 
receiver. In total, we observed 17 systems from our sample. 
The Observing log is provided in Table IA.1I For the GBT ob- 
servations of 9 Mg ii systems, we used the Spectral Processor 
as the backend with a bandwidth of 2.5 MHz split into 1024 
spectral channels. In addition, fast time sampling of 1 sec was 
used. The observations were done in standard position-switching 
mode with typically 5 min spent on-source and same amount of 
time spent at the off-source position. The data were acquired in 
two linear polarization channels i.e. XX and YY. The GMRT 
and WSRT observations were done using baseband bandwidth 
of 1MHz split into 128 and 2.5 MHz split into 1024 frequency 
channels respectively. For flux density and bandpass calibration, 
the standard flux density calibrators such as 3C48, 3C147 and 
3C286 were observed. In GMRT observations, a phase calibra- 
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0.96 



_i_ 



_l_ 



750.9 751 751.1 751.2 748.8 743.1 743.4 
J 140856 -075227 z= 1.2753 J203155+ 121941 z= 1.1 158 



0.98 - 



0.96 - 




_ 0.96 




624.4 



624.6 
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671.4 
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Fig. 2. Spectra of 21-cm absorption detections from the GBT {top) and GMRT (bottom). Individual Gaussian components and 
resultant fits to the absorption profiles are overplotted as dotted and continuous lines respectively. Residuals, with an arbitrary flux 
offset for clarity, are also shown. 



tor was also observed for lOmin after every ~45 min to get re- 
liable phase solutions. The data were acquired in the two polar- 
ization channels RR and LL. The total on-source time after ex- 
cluding telescope set-up time and calibration overheads are also 
provided in Table lA.il 



The GMRT and WSRT data were reduced using NRAO 
Astronomical Image Processing Syste m (AIPS) fo l lowin g 
the standard procedures as described in iGupta et aD (120061) . 
Unfortunately, the WSRT data for the systems towards 
J1007+2251, Jl 148+5254, J1216+5843 and J1329+1053 were 
severely affected by RFI. Although no useful spectrum could 
be obtained for these systems, the data were still useful for 
measuring the flux densities of J1007+2251, J1216+5843 and 
J1329+1053 at the redshifted 21-cm frequencies, and estimat- 
ing the accurate optical depths from the GBT data. The GMRT 
and WSRT observations resulted in 2 new 21-cm absorption de- 
tections and 21-cm optical depth limits for 6 systems. The con- 
tinuum images of two quasars that are resolved in the GMRT 
observations are presented in Fig. Q] The quasar J0956+4046 is 
resolved into two peaks. We refer to the southern component 
which is coincident with the optical QSO as PI and the other 
component as P2 in the following text. 



We used NRAO's GBTIDL package to develop a pipeline to 
automatically analyse the GBT spectral-line datasets. After ex- 
cluding the time ranges for which no useful data were obtained, 
the data were processed through this pipeline. The pipeline cali- 
brates each data record individually and flags the spectral chan- 
nels with deviations larger than 5<x seen in either XX or YY 
as affected by RFI. After subtracting a first order baseline these 
data are averaged to produce the baseline (i.e. continuum) sub- 
tracted spectra for XX and YY. The baseline fit and statistics for 
flagging are determined using the spectral region that excludes 
central 25% and last 10% channels at both ends of the spectrum. 
If necessary a first-order cubic spline was fitted to the averaged 
XX and YY spectra obtained from the pipeline, which were then 
combined to produce the stokes I spectrum. The spectrum was 
then shifted to the heliocentric frame. The multiple-epoch spec- 
tra for a system were then resampled onto the same frequency 
scale and combined to produce the final spectrum. These GBT 
observations resulted in 2 new 21-cm absorption detections and 
21-cm optical depth limits for 7 systems. 

The confusion due to other radio sources in the primary beam 
is not an obstacle in determining accurate flux densities of the 
background quasars and optical depths for the systems observed 
with the GMRT and WSRT. But for the systems observed with 
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Fig. 3. GBT, GMRT and WSRT spectra of the Mg 11 systems with 21 -cm non-detection. Arrows mark the expected positions of 
21-cm absorption lines based on metal absorption lines. For J0956+4046, PI and P2 correspond to spectra towards the southern and 
the northern peaks respectively (see Fig. [T]and the text for details). 



the GBT, the rms confusion due to other sources in the beam can 
be the dominant effect that limits the accuracy with which the 
21-cm optical depth can be determined. We used interferometric 
images from our WSRT observations and literature to address 
this issue. The details are provided in Appendix iBl and the flux 
densities and optical depth values for all the systems are listed in 
Table EU 



In summary, our spectroscopic observations have resulted in 
4 new 21-cm absorption detections that are shown in Fig. [2] The 
spectra of 13 21-cm non-detections are presented in Fig. [3] The 
GMRT spectra are at a spectral resolution of ~4kms~ 1 whereas 
GBT and WSRT spectra have been smoothed to ~4kms~ 1 . For 



the exact resolution of the spectra plotted in Figs. [2] and [3] see 
Table ED 



2.2. Details of new 21 -cm absorption detections 

Now we present the details of 4 new 21-cm absorption detec- 
tions (Fig. 0. In Table [2] we provide the details of Gaussian 
components i.e. the absorption redshift (z a bs)> me full width at 
half maximum (FWHM) and the peak optical depth (t p ) fitted to 
these spectra. These Gaussian components provide convenient 
parametrization of absorption lines even if they do not necessar- 
ily represent the actual physical components. 
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Table 2. Details of multiple Gaussian fits to 21 -cm absorption 
profiles. Listed from left to right are the names of quasars, and 
the absorption redshifts, FWHMs and peak optical depths of the 
Gaussian components. 



Quasar 


^abs 


rWrlM 


t p 






(kms- 1 ) 




J0457-2324 


0.89132 


14+2 


0.013+0.002 


J1252+4427 


0.91138 


35±6 


0.027+0.002 




0.91092 


61+16 


0.010+0.002 




0.91165 


27±6 


0.017+0.003 


J1408-0752 1 


1.27464 


11±1 


0.031+0.002 




1.27440 


9±2 


0.013+0.002 




1.27475 


10±6 


0.006+0.002 


J2031+1219 


1.11558 


78±24 


0.013+0.002 




1.11614 


62±12 


0.018+0.002 



Notes. (,) Corresponds to the spectrum taken on 21/6/2009 (see Fig.[2j- 



1 . z abs =0.8922 system towards J0457-2324: The background 
quasar, also known as PKS0454-234, is a BL Lac ob- 
j ect and highly po larized with optical polarization of 27% 
dWills et al.l 1 1992b - The 21-cm absorption detected in our 
GBT spectrum has velocity widthQ AV=23.4km s _1 and in- 
tegrated 21-cm optical depth, Jrdv=0.20+0.02km s _1 . The 
H i column density of an optically thin cloud covering the 
fraction f c of the background radio source is related to the 
integrated 21-cm optical depth and the spin temperature (T s ) 
through 

JV(Hi) = 1.823 x 10 18 — f r(v)dvcirr 2 . (1) 

fc J 

This corresponds to an H i column density of 
N(U i)=5.8xl0 19 (T s /100)(0.63// c )cirr 2 for the absorber 
towards J0457-2324. Here, we have assumed f c = Cf where 
Cf is defined as the ratio of the flux density of the VLBA 
'core' to the total arcsecond scale flux density (see Section[3] 
for the details regarding f c ). 

2. z abs =0.9112 system towards J1252+4427: The background 
quasar is resolved in FIRST survey with a deconvolved size 
of 3.68"xl.01". The 21-cm absorption detected in our GBT 
spectrum has AV=134km s'and jTdv=2.06+0.13km s _1 . 
This corresponds to N(H i)=3.8xl0 20 (T J /100)(l// c ) ctrT 2 . 
No subarcsecond scale images are available for this quasar 
to constrain f c . 

3. zabs =1.2753 system towards J1408-0755: The background 
quasar is a blazar at z- 1 .494 and the associated radio jet 
at milliarcsecond scale exhibits the highes t superluminal ve - 
locity (>25c) amongst the known blazars dPiner et al 1 l2006h . 
The 21-cm absorption detected in our GMRT spectrum has 
AV=45.0km s" 1 and JVdv=0.5 1+0.05 km s _1 . This cor- 
responds to 7Y(H i)=2.5xl0 20 (T. s /100)(0.38// c )cm- 2 , where 
we have assumed f c = Cf. We reobserved this system with 
the GMRT after 7 months (Table IA.U to obtain the higher 
resolution spectrum and investigate the variability in the 21- 
cm absorption profile due to jet proper motion. This higher 
resolution spectrum is presented in the top panel of Fig. |4] 
It shows three well detached components: Al, A2 and A3. 

2 Velocity width of absorption profile corresponding to the 5% and 
95% percentiles of the appar ent optical depth dis tribution estimated us- 
ing the method described bv lLedoux et1rilll2006h . 
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Fig. 4. 21-cm absorption spectra of Zabs= 1.2753 Mg n ab- 
sorber towards J 1408 -0752. Top panel shows the higher res- 
olution (1.9km s _1 ) spectrum taken on 22/1/2010 and bottom 
panel shows the lower resolution (3.8 km s _1 ) spectrum taken on 
21/6/2009. The 22/1/2010 spectrum smoothed to the resolution 
of 21/6/2009 spectrum is also plotted (dashed line) in the bot- 
tom panel. The shaded region corresponds to the RFI affected 
spectral range in the higher resolution spectrum. 



The detection of multiple components in this absorber could 
be related to the complex morphology exhibited by this ra- 
dio source at milliarcsecond scales. While the components 
Al and A2 are clearly detected at both the epochs, the same 
cannot be said about the component A3. The component A3 
has J rdv=0.07+0.02km s _1 and is detected with a signif- 
icance of 3.7<t at the second epoch. Unlike the features in 
the RFI affected spectral range shown in Fig. |4] the compo- 
nent A3 is also consistently reproduced in the spectra from 
the two polarisations. However, due to the lower resolution 
and S/N we cannot be certain of the presence of this compo- 
nent at the first epoch. If this component is indeed real and 
was not present at the first epoch then it would be related to 
the small scale structure in the CNM gas being tr aced by the 
jet co mponents (e.g. C2, C3 or C4 as defined bv iPiner et al.l 
20061). The proper motion of ~0.3 mas yr 1 as measured by 
Piner et all d2006l) for these components would correspond 
to a scale of ~2pc at the absorber redshift for a period of 
seven months. The follow-up GMRT observations to inves- 
tigate this further are in progress. 
4. z abs =1.1158 system towards J2031+1219: The background 
radio source is a BL Lac object at z=1.215. This radio source 
is extended in our GMRT imag e (Fig . [TJ. The VLA A- 
array map from iRector & Stockel (12001 h shows this quasar 
to be core dominated at the arcsecond scales. A broad and 
shallow 21-cm absorption is detected towards this quasar. 
The 21-cm absorption detected in our GMRT spectrum has 
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AV=143km s^'and /rdv=2.11±0.15km s _1 . This corre- 
sponds to N(H i)=8.6xl0 20 (T J /100)(0.45// c .)crrr 2 , where 
we have assumed f c = c/. 



2.3. Defining a homogeneous sample of 21 -cm absorbers: 
measurements from literature 

The purpose of this Section is to define a sample that includes 
only the 21 -cm absorption measurements of the Mg n systems 
that have been primarily selected on the basis of W r > 1 A crite- 
rion, and is therefore suitable for determining the 21-cm detec- 
tion rate in the strong Mg n systems and its redshift evolution in 
an unbiased way. Hereafter, we refer to this homogeneous sam- 
ple as S 1 . From our survey presented in this paper and G09, in to- 
tal, we have 21-cm absorption measurements for 50 strong Mg n 
systems over the redshift range 0.5 < z < 1 .5 albeit with a small 
gap at z ~ 1 • The redshift distribution of these systems is plotted 
in the Fig. [5] Only 46 of these systems that are selected primarily 
on the basis of W, > 1 A Mg n absorption should be included in 
the sample S 1 . We do not include 4 systems towards the blazars 
mentioned in Section 12.1.1 1 in SI as these were specifically se- 
lected on the basis of metal absorption line properties suggesting 
the H i column densities similar to the DLAs (see Section |5T1 >. 
It is probably not just a coincidence that 3 of these systems are 
detected in the 21-cm absorption (Section l2~!2l . 

Next we enlarge the sample SI by considering those 21-cm 
absorption measurements from the literature that match our sam- 
ple definition i.e. (1) fall in the redshift range: 0.5 < z < 1.0 
and 1.1 < z < 1.5, and (2) selected primarily on the basis 
of Wt > 1 A criterion. Since the inclusion of individually re- 
ported detections from the literature can bias our detection rate 
estimates, we tak e meas urem ents only from t he Mg n selected 
samples of lLanel <l2000h and iKanekar et all <|20Q9I) for which 
both the 21-cm detectio ns and non-detections have been system- 
atically reported. From IKanekar et al.l (120091) . the independent 
measurements for 4 systems at 1.1 < z < 1.5 and 18 systems at 
0.5 < z < 1 .0 are available. These are listed in Table IC.ll We 
have measured metal absorption line equivalent widths for these 
systems from the SDSS and the VLT archival spectra by inte- 
grating over the absorption profile in the same way as for our 
systems. Measure ments for 13 W r >1 A systems at 0.5 < z < 1.0 
are available from lLanel (12000 ) . These are listed in Table ID. fl 

Thus, we have a homogeneous sample SI of 8 1 strong Mg n 
systems with 21-cm absorption measurements. This sample hav- 
ing 12 21-cm detections will be used for investigating the red- 
shift evolution of 21-cm detection rate amongst the Mg n ab- 
sorbers. The full sample, hereafter referred as S2, of 85 sys- 
tems that also includes 4 above mentioned systems towards the 
blazars will be used when examining the relationship of 2 1 -cm 
absorption line properties with the radio structure and metal ab- 
sorption line ratios. A summary of these samples is provided in 
Table [3] In Fig. [5] we plot 21-cm optical depths as a function 
of redshift for the sample S2. For uniformity and ease of com- 
parison with the measurements from literature, the 21-cm op- 
tical depth limits ( J T3o-dv) for all the systems have been com- 
puted for a velocity resolution of 10 km s -1 . The median J T3o-dv 
achieved for the systems from [Lane (2000) is ~0.64km s _1 . 
This is at least a factor three lower sensitivity as compared to 
that achieved for ou r measurements or for the systems from 
IKanekar et ail d2009l) . The method we use to estimate the 21- 
cm detection rates by taking into account the variation in opti- 




Fig. 5. Redshift distribution of 21-cm optical depths for the 
strong Mg n absorb ers from this paper a nd G 09 (circles), 
IKanekar etail d2009l) (squares) and lLanel (120001) (triangles). 
The filled symbols are for 21-cm absorption detections and the 
open symbols correspond to 3<x optical depth limits for 21-cm 
non-detections. In the top panel, dotted, dashed and solid line 
histograms are for the sample S2 with J T3o-dv<0.1, 0.3 and 
0.5 km s respectively. 



cal depth sensitivities towards different quasar sight lines is de- 
scribed in Section |3~T1 



2.4. VLBA continuum observations 
2.4.1. Motivation and sample 

The detectability of 21-cm absorption, in addition to the optical 
depth, depends directly on the covering factor, f c , of the absorb- 
ing gas (see Equation Q]). Since radio sources in general exhibit 
emission over a wide range of physical scales ranging from pc 
to kpc the f c will depend both on the overall extent of the back- 
ground quasar and the small-scale structure (i.e. dumpiness) of 
absorbing gas. The motivation behind these VLBA observations 
is to understand the relationship between the detectability of 21- 
cm absorption and the structure in the radio source and absorbing 
gas at the scales of tens of pc. For this purpose we selected all 
the systems from our sample that are towards the quasars com- 
pact at arcsecond scales. This list of 42 Mg n systems towards 
40 quasars is drawn from Table 1 of G09 and Table Q~|of this pa- 
per, and it includes all the systems towards the compact quasars 
from our sample with 21-cm absorption spectra. Exceptions are 
J0457-2324 and J2031 + 1219. The former, although compact at 
arcsecond scales, could not be observed whereas observations 
of BL Lac J203 1 + 1219, extended at arcsecond scales, could be 
included within our observing schedule. 
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Table 3. Summary of all the W r > 1 A Mg n absorbers at 0.5 < z < 1 and 1 . 1 < z < 1 .5 with 2 1 -cm absorption measurements considered 
here. The samples SI and S2 discussed in the text are comprised of the subsamples B-E and A-E, respectively. 



Sample 


Number 


21 -cm absorption 


Selection 


Reference 




of systems 


detections 


criterion 




A 


4 


3 


W s >1A + Zn n/Crii/Mnii 


This paper 


B 


13 


1 


Wr >1A 




C 


32+ 


8 




GuDta et al. (2009) 


D 


23 


1 


»> 


Kanekar et al. (2009) 


E 


13 


2 




Lane (20001 



Notes. (,) Measurements of 33 W T >1 AMg n systems are available from G09. However, the optical depth value for the system towards Jl 145+0455 
is taken from the more sensitive spectrum of iKanekar et alj d2009h . 



2.4.2. Observations and data reduction 

The VLBA 20-cm observations were done as part of a larger sur- 
vey to obtain milliarcsecond scale images of the quasars with 
DLAs and Mg n systems along the sight line. Observations 
were carried out on 28/10/2008 (~5hrs), 7/11/ 2008 (~14hrs) 
21/2/ 2010 (-llhrs) and 10/6/2010 (~17hrs) ( ISrianand et all 
l2012h . We used eight 8 MHz baseband channels i.e. the total 
bandwidth of 64 MHz. Each baseband channel was split into 32 
spectral points. Both the right-hand and left-hand circular polar- 
ization channels were recorded. Two bit sampling and the time 
resolution of 2 seconds were used. 

The observations were done using nodding-style phase- 
referencing with a cycle time of ~5 min i.e.~3 min on the source 
and ~ 1.5 min on the phase-referencing calibrator. The phase- 
referencing calibrators were selected from the VLBA Calibrator 
Survey (VCSfl at 13-cm and 3.6-cm. In order to improve the 
uv-coverage, the total time on each source was split into snap- 
shots over a number of different hour angles. Each source was 
typically observed for a total of ~30min. Strong fringe find- 
ers/ bandpass calibrators such as J0555+3948, J0927+3902, 
J1800+3848 and J2253+1608 were also observed every ~3hrs 
for 4-5 min. 

The data were calibrated and imaged using AIPS and 
DIFMAP in a standard way (e.g. lBeaslev et alj2002l) . The global 
fringe fitting was performed on the phase-referencing calibra- 
tors. The delays, rates and phases estimated from these were 
transferred to the target source which was then self-calibrated 
until the final image was obtained. The VLBA maps of these 40 
quasars are presented in Fig. IE. II 

The milliarcsecond scale structure was characterised by fit- 
ting Gaussian models to the self-calibrated visibilities. The re- 
sults of the model fitting are provided in columns#5-ll of 
Table [4] We define /vlba as the ratio of total flux density de- 
tected in the VLBA and the arcsecond scale image. Similarly, 
Cf is defined as the ratio of the flux density of the VLBA 'core' 
component to the total arcsecond scale flux (see Section|3]for de- 
tails). The largest linear size (LS) represents separation between 
the farthest components of the radio source at the absorber red- 
shift. For the compact sources represented by single component 
we take the major axis of the deconvolved component as the up- 
per limit on source size. The /vlba and cj are discussed in detail 
in the following sections where we use them to estimate f c and 
correct observed 21 -cm optical depths for partial coverage. 



2.5. Milliarcsecond scale images from literature and covering 
factor 

The radio morphology of all the quasars in the 21 -cm absorption 
sample S2 are important to our later discussion. Therefore, for 
the quasars that were not covered in our VLBA 20 cm survey we 
have compiled this information from literature. For J0457 -2324, 
Jl 145+0455 and J2129-1538 from our sample, the c f , /vlba 
and largest linear size (LS) obtained using maps at 13-cm from 
the VCS are provided at the end of Table |4] 

In addition, the milliarcsecond imag es for the 14 q uasars 
with 15 intervening systems taken from IKanekar et all (120091) 
are available at 13-cm from the VCS, and for another quasar 
B0812+332 from our VLBA survey. The /vlba, Cf and LS for 
these sources are given in the last three columns of Table IC.fl 

From the sample of 13 systems taken from lLanel (120001) . we 
will consider here only six systems that have been observed with 
the optical depth sensitivity, Jr3o-dv<0.6km s _1 . From this, the 
milliarcsecond images at 13-cm are available for B0235 + 164, 
B04 54+039 and B 12 18+3 39 from VCS, and for B1629+139 
from Dallacasa et al ] (ll998h (TablelrTTTi. As we shall see later in 
the Section l4~Tl the remaining 7 systems that are all 21 -cm non- 
detections, have been observed with the optical depth sensitivity 
that is too poor to address the covering factor related issues. 

In total, we have information on milliarcsecond scale radio 
structure for 54 quasars with 57 intervening systems from the 
sample S2. The VLBA information for the 34 quasars (with 36 
systems along their line of sight) are based on 20-cm VLBA im- 
ages and for the remaining quasars it is based on 13-cm VLBA 
images. Unlike using the source sizes and flux densities obtained 
from the Gaussian parametrization of 20-cm maps from our sur- 
vey, the fvLBA and c/ at 13-cm have been estimated using the 
total and peak flux densities from the VCS. However, these val- 
ues of the Cf and /vlba for the cases with the 'core' identifica- 
tions, especially in the 13/18 sightlines where the background 
quasar can be represented by a single component, should not 
diffeiQ by more than ~10% relative to the estimates from the 
Gaussian parametrization technique. The total arcsecond scale 
flux densities at 13-cm in all these cases have been estimated 
by interpolating flux densities available from the NASA extra- 
galactic database. The upper limits on the size of quasars that 
are represented by single component have been estimated using 
the major axis of point spread function. These are conservative 
compared to those obtained using the Gaussian deconvolution. 



http://www.vlba.nrao.edu/astro/calib/ 



Based on the data from our 20-cm VLBA survey. 
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Fig. 6. Spectral index (a) vs linear size (in pc) at the absorber 
redshift for the quasars with the milliarcsecond scale images. 
The a are from the arcsecond scale images. Systems towards 
the quasars with no 'core' identifications are plotted as squares. 
The point for the system towards the quasar B 1 629 + 1 20 that has 
a=+0.54 and LLS=7.2 kpc is not plotted to keep the x-axis short. 
The histogram for the a is also shown. The filled symbols and 
shaded histogram are for 21 -cm detections. 



3. Estimating the gas covering factor and 21 -cm 
absorption detection rate 

Radio source structure often consists of a combination of com- 
pact components and extended diffuse emission. The latter due 
to its low surface brightness sensitivity will only contribute to 
the 2 1 -cm optical depth sensitivity if the large scale structure of 
the absorbing gas is well aligned to the radio emission of the 
quasar. Such fortuitous combinations are rare. Therefore by us- 
ing the flux densities from our spectroscopic observations that at 
the best have a spatial resolution of ~30-40 kpc at the absorber 
redshift, in general, we are overestimating optical depth sensi- 
tivities. 

To show that this is indeed the case with the measurements 
considered in Table [3] in Fig. [6] we plot the spectral index0 (a) 
vs linear size (LS) for all the quasars from the sample S2 with 
milliarcsecond scale images. These spectral indices have been 
estimated between 20-cm and 50-cm/90-cm using the flux den- 
sities from the maps with arcsecond scale resolution and range 
from -1.2 to 0.7. It is clear from Fig.[6]that the 21 -cm detection 
rate is higher towards the quasars with the flat/inverted spectral 
indices. The 44^% (7/16) of systems towards the quasars with 
a <0 from the sample S2 are detected in 21 -cm absorption as 
compared to the 15 6 % (6/41) towards the quasars with a >0. 
It is also clear that 73% of the 21 -cm detections are towards 
the quasars with LS<100pc and amongst the quasars with LS < 
lOOpc the detection rate is less towards the quasars with steeper 
spectrum (i.e a > 0). This may imply that the characteristic size 
of the absorbing gas is less than lOOpc. The low detection rate 



We adopt the definition that the flux density at frequency v, S v °cv °. 



in the quasars with the steeper spectrum will then be related to 
the poor optical depth sensitivity to detect 21 -cm absorption to- 
wards the 'cores'. The term 'core' is used here to represent the 
radio component coincident with the optical QSO. 

Thus, identifying the 'core' component in the sample S2 is 
important to confirm these suggestions and estimate f c of the 
absorbing gas. For this, we need milliarcsecond scale maps at 
multiple frequencies to estimate the spectral indices of the var- 
ious components detected in the 20-cm or 13-cm VLBA maps 
discussed in Section 12.51 For the 42 systems, images of the 
background quasar at more than one frequency are available 
through our observations, the VCS and the VLBA Imaging and 
Polarization Survey (VIPS). Using these images at 20-cm, 13- 
cm, 6-cm and 3.6-cm we find that in 36 cases the dominant 
VLBA component has a flat spectrum at milliarcsecond scales 
and can be identified as the 'core'. Although for the quasars 
J0154-0007, J0845+4257 and J0852+3435 only images at 20- 
cm are available but all three have more than ~70% of the arc- 
second scale flux contained in the dominant VLBA component. 
We identify this component as the 'core' on the basis of a<0 
estimated at the arcsecond scale. Therefore, in total, for the 39 
systems we identify the dominant VLBA component of the back- 
ground quasar as 'core'. 

In Fig.|6l we plot 17 quasars with 18 intervening systems that 
have no 'core' identification as squares. In 6 cases, the multifre- 
quency VLBA images are available, but the dominant compo- 
nent at the milliarcsecond scales has a steep spectrum and can- 
not be identified as 'core'. No multifrequency VLBA images are 
available for the remaining 11 quasars with 12 intervening sys- 
tems. These 17 quasars with no 'core' identifications also have 
on the average steeper spectral index. The cy values for these are 
marked as '*' in Tables l4l and IC.ll This subset of 18 systems 
includes three 21-cm detections which are towards J1058+4939, 
Jl 145+0455 and J1623+0718. The former two exhibit complex 
morphology at milliarcsecond scales. The third is compact in 
20-cm VLBA map but unlike most of the other quasars that are 
compact at milliarcsecond scales has a low cy value of 0.42. We 
will discuss the nature of these systems in detail in the subse- 
quent sections. 

Thus, we have 'core' identification for 39 quasars from the 
sample S2. We will consider the following two scenarios and 
use Cf and /vlba to estimate f c for the 21-cm absorption sys- 
tems towards these. In the first scenario we define f c — cj. The 
assumption in this case is that the 21-cm absorbing gas only cov- 
ers the 'core'. If the filling factor of the absorbing gas at these 
scales is unity, this provides a conservative lower limit on the 
actual value of f c . In the second scenario we use f c = /vlba - The 
assumptions in this case are that the 21-cm absorbing gas ex- 
tends well beyond ~300pc to cover all the milliarcsecond scale 
components and the diffuse extended (>kpc scale) emission that 
got resolved out in the VLBA maps does not contribute to the 
optical depth sensitivity (see also Section l6~3l l. 

3.1. Detection rate 

In this Section, we determine 21-cm detection rates in the strong 
Mg n absorbers. Due to the reasons mentioned in Section 12.31 
we will use the homogeneous sample S 1 (Table [3} for this pur- 
pose. In the sample SI, there are 40 systems at 0.5 < z < 1 
and 41 at 1.1 <z< 1.5. There are 3 and 9 detections in the 
low- and high-z bins respectively. The top panel of Fig.[5]shows 
the distribution of optical depth sensitivities in these two redshift 
bins. We will determine the 21-cm detection rates following the 
method described in Section 7 of G09 that takes into account the 
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variation of optical depth sensitivities towards different quasar 
sight lines. The detection rate for 21-cm absorption, C, (defined 
as the fraction of Mg n systems that show detectable 21-cm ab- 
sorption with the integrated optical depths T21 greater than the 
3cr limiting value To) for the 0.5 < z < 1 and 1.1 < z < 1-5 
bins are given in Table [5] Since the lower limiting values of To 
select the weaker 21-cm absorbers and higher values select the 
stronger absorbers, in Table|5] we give detection rates for a range 
of To values suitable for the sample considered here (see the his- 
tograms in Fig.O. 

It is clear from Table [5] that the 2 1 -cm detection rate C for 
any value of T" is only slightly lower for the low-z bin. Since 
the number of detections involved are small in both the redshift 
bins, the errors are large and the noted difference is not statisti- 
cally significant (i.e say at more than 3cr level). The maximum 
difference is seen for T" o =0.1 km s _1 and is only significant at 
the level of ~0.6cr. To test this further, we perform the general- 
ized rank correlati on test that takes into account the upper limits 
dlsobe et al.lll986l) . We find a correlation coefficient of 1.37 and 
the probability that this correlation occurs by chance of 17%. 
Therefore the apparent trend is also hinted by the correlation test 
but once again without high statistical significance. 

Thus we conclude that the 21-cm absorption detection rate 
'C for the strong Mg 11 absorbers is constant over 0.5< z <1.5 
that is ~30% of the age of Universe. As the radio emission from 
quasars originate from a much larger spatial extent than the op- 
tical emission, one must pay attention to the issues related to 
radio source structure. Therefore, before discussing the evolu- 
tion of 'C any further, we investigate the impact of the radio 
structure of background quasar on the 21-cm detection rate. 




0.2 0.4 
/ T 3„ dv (kms- 1 ) 



Fig. 7. Distribution of 21-cm optical depth sensitivity for the 
quasars that are compact and extended at arcsecond scales in the 
sample S2. Shaded histogram shows the sensitivity achie ved in 
case of 21-cm absorption detections. The 7 systems from iLanel 
(I2000h with J T3o-dv>0.6 km s 1 are placed in the last bin. 



sumptions of the gas filling factor and the extent of absorbing 
gas (cf. Section[3]i. 



4. 21-cm absorption and radio source structure 

ICurranl d2010h has used angular diameter distance ratios 
(DA fl /„/DA 9 . TO ) at the absorber and background quasar redshift 
to investigate the effect of line-of-sight geometry on the cover- 
age of background radio source by the 21-cm absorbing gas. He 
suggests that since the background radio source can be more ef- 
fectively covered at the smaller DA a b s /DA CJSO , the larger values 
(>0.8) of DA a i, s /DA qso at z a bs^l can explain the lower 21-cm 
detection rates at the higher redshifts. 

In the homogeneous sample SI, only 4 out of 81 systems 
have 0.6 <DA fl /„/DA 9 . so < 0.8 and two of these are detected 
in 21-cm absorption. The remaining 95% of the Mg 11 systems 
have DA a £ )J /DA, ?JO >0.8 and dominate the statistics. As expected, 
all the 4 systems with DA a b s fDA qso < 0.8 are at z a bs<l- For 
7~o=0.1 km s , they correspond to C=0.50^°. Interestingly, this 
is same as the C for 7~~o=0.1km s _1 we obtained in the pre- 
vious Section for the z > 1 bin with all the systems having 
DA fl /„/DA 9J0 > 0.8. Therefore, we do not see an y significant in- 
fiuence of the geometric effects as quantified bv lCurranl (1201 Ol) 
on the detectability of 21-cm absorption in our sample. This also 
implies that the redshift evolution of 21-cm detection rates over 
the low- and high-z bins in Table [5]is not affected by any evolu- 
tion in DA abs fDAg S0 . 

However, this does not mean that the 21-cm detection rate 
estimates are not affected by the partial coverage of the back- 
ground quasar. As discussed in Section [3] the 21-cm optical 
depth sensitivities from the GBT, GMRT and WSRT observa- 
tions are overestimated. Our purpose here is to estimate the 21- 
cm optical depth sensitivities and detection rates corrected for 
partial coverage. For this we will use f c as estimated from the 
arcsecond and milliarcsecond scale images under reasonable as- 



4.1. Arcsecond scale radio structure 

The sample S2 has 56 and 31 systems towards the quasars that 
are compact (deconvolved sizes <2") and extended at arcsecond 
scales, respectively. The distribution of optical depth sensitivities 
for these are plotted in Fig. [7] The 21-cm absorption detection 
rates for T" o =0.3km s _1 towards the compact and the extended 
quasars are 10^% and 14^' 3 % respectively. The detection rates 
are also similar for T=0.1km s and 0.5 km s _1 , and do not 
show any significant difference when the sample is split into the 
two redshift bins as given in Table [5] However, it is interesting 
to note from Fig. Qthat while in the case of systems towards the 
extended quasars the 21-cm absorption is detected only when 
f T3o-dv< 0.2 km s -1 , in case of compact quasars the 21-cm ab- 
sorption is also detected for J T3o-dv>0.4 km s _1 (without any 
correction for the partial coverage). To understand this, we now 
discuss these detections in detail. 

Amongst the systems towards the compact quasars, the three 
21-cm detections in the tail of J T3o-dv distribution are towards 
J0850+5159, J0852+3435 and J2340-0053. The optical depth 
sensitivities achieved are lower as these systems are towards rel- 
atively weaker radio sources in our sample. Still, the 21-cm ab- 
sorption is detected due to the large 21-cm optical depth. The 
rare 2 175 A UV-bump is detected in the first two cases suggest- 
ing the presence of cold gas with a large c olumn density of dust 
(ISrianand et al.l2008tlKulkarni et alj|201 ll see also SectionEHli. 
In the third case, the 21-cm line width is consistent with the ab- 
sorbing gas having T s <200K resulting in large 21-cm optical 
depth (G09; Rahmani et al. in prep.). 

Now we discuss the 21-cm detections towards the quasars 
that are extended at arcsecond scales. The 5 detections in 
this case are towards B0235 + 164, J0804+3012, Jl 145+0455, 
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Table 5. The fraction of 21-cm absorption detections (C) and the number density of 21-cm absorbers («2i) for the different 21-cm 
optical depth limits (To) using the homogeneous sample SI. N is the total number of systems in the sample with J r^wdv < To, 
and N21 is the number of 21-cm detections with f T^iodv < To and J rdv > To- 



z range 


Sample 


T 
(km s _1 ) 


Systems (N) 


Detections (N21) 


C 


"21 


1.1-1.5 


Ours (i.e. This paper + G09) 


U. 1U 






4 


n cn+0.40 
0.50_ 024 


n 1 •3+0.11 






0.20 


22 


4 


n 1 o+O 14 


n nc+0 04 

°- 05 -RR? 






0.30 


27 


3 


0.1 1 


03 + ™ 






0.50 


36 


4 


r\ 1 1 4-0 fit) 

0. 1 1 


°-°3-o:m 


1.1-1.5 


Ours + literature 


0.10 


10 


5 


0.50JI 


^» y ^ _i_n no 

13+u.iw 

-R-flG 






0.20 


26 


5 


n 1 n+0 13 


n nc+0 04 
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0.30 


32 


4 


13+™ 
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41 


5 


1 2 +u,us 
-RQA 


03 

"■"-'-o.oi 


1.1-1.5 


Ours + literature with VLBA 


0.10 


4 


2 


50 +i "" 


13 +017 




(assuming f c = c f ) 


0.20 


8 


1 


r\ 1 4-0 9Q 

U,AJ -Q.10 


03 +0 07 

"•"■'-0.03 






0.30 


10 


1 


10 


03 +0 06 
"'"•'-002 






0.50 


16 


2 


n 1 1+6 17 

1 1 > TV. 1 / 

u - -0.08 


f\ r\o +fi 04 

a° 3 -o:o2 


1.1-1.5 


Ours + literature with VLBA 


0.10 


4 


2 


so +u ; '" 

-0.32 


n -t s% _i_n 1 1 

1 3+" 1 ' 
-Q-Q9 

03 

-Q 03 




(assuming f c = f VLBA ) 


0.20 


8 


1 


1 

-Rlfl 






0.30 


13 


1 


0.08!»; 


„ „ _ .A its 

02 

-R-R2 






0.50 


17 


3 


<-v 1 n4-0 18 


r\ r\/--i-Pl OS 

05 

u - UJ -0.03 


0.5-1.0 


Ours 


0.10 


4 


1 


r, c-4-0 18 

() 2v 

-Q-2JL 
n 1 -7+0 38 

°-ul 


n n^+n in 
04 

-Q-Q4 






0.20 


6 


1 


n m+fi 07 

0.03!""' 






0.30 


8 


1 


n 1 o+O 29 

°- 13 -R42 


0.02!""' 

02 

03 + ° fe 
-Q-Q2 
n no +002 

0.02_ 00| 






0.50 


9 


1 


lr"- 2 ' 5 

R-9§ 


0.5-1.0 


Ours + literature 


0.10 


20 


3 


is+ - 15 

-Q-OS 






U.zU 


23 


-> 

i 


n 1 o+O 12 






0.30 


28 


3 


11+ 010 


02 +002 

-8R 






0.50 


33 


2 


o.o6:| 


O01+001 

-Q01 


0.5-1.0 


Ours + literature with VLBA 


0.10 


5 


2 


o.4o:»| 


07 +m 
-SR6 




(assuming f c = c/) 


0.20 


10 


2 


020+025 


04 + ° "' 

-QQ3 




0.30 


12 


2 




'-0.14 


03 +i5 " 4 






0.50 


15 


2 


13 +a2 ° 
-Qll 


02 +0 03 


0.5-1.0 


Ours + literature with VLBA 


0.10 


5 


2 


0.40!;* 


007+0.10 

u - u '-0.06 




(assuming f c = f VLBA ) 


0.20 


1 1 


2 


0.18!; 
0.14 +iii > 


03 +0 05 

"•"-'-0.03 






0.30 


14 


2 


03 +0 04 

-QQ2 

02 +m 
"■"--00? 






0.50 


15 


2 


13 +8:1^ 



J1252+4427 and J2031 + 1219. The subarcsecond scale images 
are available for 4 of these quasars (i.e. except J1252+4427) that 
are also amongst the brightest radio sources in our sample. As 
expected the estimated correction to account for the partial cov- 
erage is also large (cy <0.5) for these. Despite this the 21-cm 
optical depth sensitivity achieved towards the 'core' component 
is exceptionally good (i.e. Jr3o-dv<0.2 km s~'even after correct- 
ing for the f c assuming f c = cf) for these 4 systems. The good 
optical depth sensitivity towards the milliarcsecond scale 'core' 
coupled with the large f c correction (cj <0.5) can be taken to 
imply that these quasars also have good optical depth sensitivity 
at the physical scales beyond (>30pcf3 the 'core' component. 
We notice that the 21-cm absorption in all these 5 cases exhibits 
large velocity widths i.e. AV=100-150km s _1 . Compared to this 
the widths of absorption lines towards the compact quasars plot- 
ted in the top panel of Fig. Q are less than 100km s~'. This 
favours the idea that larger AV observed towards the extended 
sources arise from the 21-cm absorbing gas that has transverse 
extent of >30pc. We examine this further in the next Section. 



6 Based on the typical resolution at the absorber redshift achieved in 
the VLBA images used for 'core' identification. 



4.2. Correlations with VLBA morphological parameters 

Now we explore the correlations between the VLBA morpho- 
logical and 21-cm absorption parameters using the sample S2. 
In Fig. [8]we plot J rdv/fvLBA as the function of linear size. If we 
consider only those systems that are towards the quasars with 
the clear 'core' identification (plotted as circles in Fig. |5J, we 
see that almost all the detections (except one) are towards the 
quasars with the linear size, LS<100pc. However, the gener- 
alised rank correlation test suggests no correlation between the 
linear size and the Jrdv/c/, or J rdv/fvLBA- The correlation be- 
tween the linear size and J rdv/fvLBA is only significant at the 
level of 0.4cr with a probability that it can arise due to a chance 
of 7 1 % . The significance is even lower (0. 1 cf) for the correlation 
between the linear size and J rdv/c /. 

Now we return to the issue of relationship between the radio 
structure at subarcsecond scales and the width of 21-cm absorp- 
tion lines. In Fig. [9] we plot the width of 21-cm absorption lines 
as the function of linear size measured from the VLBA maps 
for all the 21-cm detections from the sample S2. The only ex- 
ception is the system towards J0804+3012 for which the lin- 
ear size is measur ed from the subarcsecond scale images of 
iKunert et al.l d2002l) . There are four systems in this figure show- 
ing AV>100km s~'. These four are identified in Section |4~T1 as 
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log Linear size (pc) 



Fig. 8. Integrated 21 -cm optical depth corrected for partial cov- 
erage assuming f c = fvLBA vs linear size (in pc) at the absorber 
redshift for the sample S2. Systems towards the quasars with 
no 'core' identifications are plotted as squares. 21 -cm detec- 
tions are shown as filled symbols where 21 -cm non-detections 
are represented by open symbols. The point for the system to- 
wards B1629+120 with LS=7.2kpc is not plotted to keep LS- 
axis short. 
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Fig. 9. Width (AV) of 21-cm absorption line in km s _1 is plotted 
against the linear size in pc measured from the VLBA maps. 
The systems towards the quasars with the 'core' identification 
are plotted as circles and the remaining are shown as squares. 



the ones towards the radio sources having extended morphology 
at arcsecond scales. It is also interesting to note that there are 
two systems with AV<40km s _1 towards the quasars that show 
extended morphology in the VLBA images (i.e J1058+4939 and 
B1629+120). The correlation analysis suggests only a weak cor- 
relation at the level of 1.54cr with the probability that it can be 
by chance of 0.12. 

The lack of strong correlation between AV and projected lin- 
ear size is not surprising. This is because, even if we only look 
at compact radio sources we expect a spread in the line of sight 
AV as measured by the distribution of Mg n equivalent widths. 
Such line of sight velocity spread will dilute any possible cor- 
relation between AV and LS that is expected to originate from 
the velocity gradients in the transverse direction as traced by the 
extended radio sources. In our sample, we do not find any cor- 
relation between W r and AV. The correlation analysis suggests 
a significance of 0.27<x with a chance probability of 0.74. Note 
that Curran (2010) also finds only 1.73<x correlation between W r 
and AV for a sample comprising both DLAs and W r >1 A i.e. 
strong Mg n absorbers. 

In summary, we find a large fraction of 21-cm detections to 
be towards the quasars having linear size < lOOpc in the VLBA 
images. The largest velocity widths are seen towards the quasars 
that show extended structure at arcsecond scales. However, we 
do not find strong correlation between either J rdv or AV with 
the linear size measured from the milliarcsecond scale images. 
All this can be understood if the absorbing gas is patchy with a 
typical correlation length of ~30-100pc. It is quite likely that the 
narrow absorption components with the larger optical depths are 
associated with even smaller clouds (Srianand et al. in prep.). 



4.3. Milliarcsecond scale structure and 21 -cm detection rate 

In this Section, we determine the 21-cm detection rates for the 
strong Mg n absorbers after correcting the 21-cm optical depths 
for partial coverage. For this purpose, we will use the sample S 1 
that only consists of the systems primarily selected on the basis 
of W, > 1 A criterion as described in the Section 12.31 As noted 
above, if the 21-cm absorbing gas indeed extends beyond 30 pc 
then under the assumption of unit filling factor we can use c j and 
fviBA to estimate f c and obtain a realistic estimate of the optical 
depth sensitivity. As discussed in Section. [3] we have two possi- 
ble values of f c for 35 systems with 'core' identifications in the 
sample SI. Just to recollect, in the first case we assume f c = c/. 
The f c in this case provides a conservative lower limit on the 
actual value of f c . In the second case of f c - /vlba, we assume 
that the absorbing gas extends over 300 pc to cover all the VLBA 
components. We determine the 21-cm absorption detection rates 
for these two scenarios using the corrected optical depth sensi- 
tivity limits. These detection rates for the different values of To 
are summarised in Table 

It is clear from Table [5] that the detection rates in the low-z 
bin are slightly higher in this case compared to the results ob- 
tained with no correction for the partial coverage. The slight gap 
noted between the low- and high-z detection rates, when no cor- 
rection for partial coverage is applied, is also reduced. The max- 
imum increase in C after f c correction is by a factor 2.7 for the 
low-z bin. Thus covering factor is an important issue in measur- 
ing the redshift distribution of 21-cm absorption systems accu- 
rately. One may underestimate detection rate by up to factor ~2 
due to this effect alone. 
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Fig. 10. Ri and R2 vs DR are plotted for the systems towards quasars with the VLBA images. Dashed lines in the panels are for 
median values of DR=1.1, Ri = 1.47 and R2=0.27. Systems at 0.5 < z < 1 are plotted as squares and 1.1 < z < 1.5 as circles. Filled 
symbols are for 21 -cm de tections. Values for system towards J0808+4950 (G09), J1058+4939 (G09), Jl 145+0455 (Table|OJ and 
B1629+120 (lLan el l2000h are labelled. Median error bars are shown at the top-left corner. System towards J2129-1538 is omitted 
due to large uncertainties on equivalent widths. 
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Fig. 11. Distribution of the number of Mg 11 systems for different subsamples from the sample S2 with the 'core' identification. 
Systems at 1.1 < z < 1.5 are plotted in the top panel and systems at 0.5 < z < 1 in the bottom panel. Shaded histogram corresponds 
to detections and optical depths have been corrected for partial coverage assuming f c = c/. For each subsample four bars correspond 
to no-cutoff, T = 0.3, 0.2 and 0.1 km s~' respectively going from left to right. Detection rate, C, in each case is given at the top of 
bar. 
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5. 21 -cm absorption and properties from their 
optical spectra 

5.1. Metal absorption lines 

It has been well established that the fracti on of Mg n sys - 
tems that are DLAs increases with the W x dRao et alj l2006l) . 
However, there exists a large scatter in the W r and N(H i) rela- 
tionship, and W r cannot be used as an indicator of the N(H i) 
dRao et alj |2011|) . On the other hand, the Mg n doublet ra- 
tio (DR), W(Mg ii/12796)/W(Fe H/12600) (defined as RO and 
W(Mg Ll2852)/W(Mg H/12796) (defined as R 2 ) can be used as 
a robust in dicator of t he pre sence of high N(H i). From the HST 
sample of iRao et all d2006h . we notice that 27^% and 38^% 
of the strong Mg n systems are DLAs at 0.5 < z <1.0 and 
1.0< z <1.5 respectively. The DLA fraction can be as high as 
50% if further restrictions on the DR, Rj and R2 are imposed. 
G09 also noticed that a higher 21 -cm detection rate can be ob- 
tained with the systems limited to a restricted range of these pa- 
rameters. 

In order to confirm if the 21 -cm absorption detections oc- 
cupy any preferred location in the parameter space defined by 
the metal absorption line ratios, we plot in Fig.[l0]the ratios Ri 
and R2 as a function of DR. The Mg 1, Mg 11 and Fe n equiv- 
alent widths for the systems taken from iKanekar et al.l d2009h 
(see Table [CTTb have been estimated using the SDSS and VLT 
archival spectra following the same procedure as done in G09 
and TableQ] The equivalent widths for the systems that are from 
Lane (l2000t) are given in Table ID. fl These have been taken from 
Land d2000t) . When comparing the properties of the radio and 
optical absorption lines it is important to be able to address the 
issues that could arise due to the different gas volumes being 
probed by the radio and optical sight lines. For this reason, from 
the sample S2, we will consider here only the Mg 11 systems to- 
wards the quasars with the VLBA images. 

It is clear from Fig. [10] that if we split the sample along 
the median values of DR, Ri and R2, plotted as dashed lines, 
then -75% of the 21-cm detections have DR<1.1, Ri <1.47 and 
R2 >0.27. The detection rate amongst the Mg 11 systems in the 
1.1 < z < 1.5 sample, without any optical depth cut-off, is as 
high as 80% for these values of DR, R; and R2. This trend for 
1.1 < z < 1.5 was previously noted by G09. Since 75% of the 
Mg 11 systems in the sample S2 at 1.1 < z < 1.5 considered 
here come from G09, this finding is not surprising. This result 
implies that the probability of detecting cold gas via 21-cm ab- 
sorption is higher in systems with the high values of A^(H 1), 
and is consistent wit h the 3<x leve l corre lation between N(H 1) 
and T s /fc noted by ICurran etal ] (120101) . A similar trend has 
also been noted f or 21-cm detect i ons in a sample of 28 bonafide 
DLAs at z > 2 bv lSrianand et al.l d2012l) . 

In this context, it is interesting to recall that in the 
Section |2l we also presented 21-cm absorption measurements 
of the 4 Mg 11 systems i.e. towards the blazars J0457-2324, 
J1408-0752, J2031 + 1219 and J2129-1538, that were specifi- 
cally selected on the basis of the metal absorption line prop- 
erties that suggest that these systems may have H 1 column 
densities similar to the DLAs. The first three not only have 
W(Mg ii/i2796)/W(Fe ii/12600)<3/2, a criterion satisfied by the 
DLAs, they also have strong Mg 1 and Mn 11 absorption. The 
first two systems also show absorption lines of Zn 11 and Cr 11. 
In the case of system towards J2129-1538, the absorption lines 
of Na 1/1/13303.3,3303.9 are also detected. Detection of these 
species together with the Mg 1, Mg 11 and Fe 11 suggest that the H 1 
column densities should be significantly higher than the thresh- 



old used to define DLAs. Therefore, high 21-cm detection rates 
towards these 4 blazars could be due to these systems having 
high N(U 1). 

Now we discuss the two deviations from the above men- 
tioned trends. First, there are 3 21-cm detections in Fig. [TUJ la- 
belled by the quasar name, that do not satisfy both the DR-Rj 
and DR-R2 values satisfied by the other 21-cm detections. The 
doublet ratio for system towards J1058+4939 is only 1.5<x from 
the median dashed line of DR=1.1. It is interesting to note that 
the quasars J1058+4939, Jl 145+0455 and B1629+120 exhibit 
multiple radio components at sub-arcsecond scales. The dom- 
inant VLBA component for these three quasars is not 'core'. 
Therefore it is possible that the radio and optical sight lines in 
these cases are probing different gas volumes resulting in the de- 
viations seen in Fig.fTUl 

Second deviation is that at 0.5 < z < 1 only 20-30% of 
the systems with the DR< 1.1, R] <1.47 and R 2 >0.27 are de- 
tected in the 21-cm absorption whereas at 1.1 < z < 1.5 this 
fraction can be as high as 80%. This cannot be attributed to the 
lower optical depth sensitivities achieved in the low-z bin or the 
evolution of H 1 column density amongst the strong Mg 11 ab- 
sorbers as the fract ion of strong Mg n systems that are DLAs in 
the HST sample of lRao et all d2006l) are same at 0.5 < z < 1 and 
1.1 < z < 1.5. To examine this in detail we plot in Fig.Qj] the 
distribution of systems in the sample S2 for the different subsets 
of Mg 11 systems defined on the basis of DR, Ri and R2. For 
each subsample, the four histograms provide the number of sys- 
tems and 21-cm detections for no optical depth cut-off, and T 
= 0.3, 0.2 and 0.1 km s _1 respectively. We have considered here 
only the subset of 39 systems with clear 'core' identifications, 
and the optical depths have been corrected for partial coverage 
assuming f c = Cf. The trends observed in Fig. [TTJare similar if 
we consider all the 56 systems with the VLBA maps and assume 
fc = fvLBA to correct for the partial coverage. 

We stress here that the metal absorption line ratios measured 
from low resolution spectra such as SDSS considered here pro- 
vide the properties averaged over several absorption components 
for a system. 21-cm absorptio n often arise from certain specific 
components (see e.g. Fig. 8 of lGupta et alll2009l) . Also the error 
bars on individual DR, Ri and R2 measurements are large (see 
Fig. [Toll, and the observed deviations could simply be due to sta- 
tistical fluctuations. Higher resolution optical spectra and most 
importantly N(H 1) measurements are required to address these 



5.2. Dust content 

Most of the Mg 11 systems and DLAs in SDSS trace diffuse ISM 
and rare ly show reddening due to dust in the spectra of individual 
quasars dBudzvnski & HewettfeOllllWild et al.l2006l:lYork et al.l 
l2006h . However in the 2 1 -cm absorption sample discussed in this 
paper, there are 4 Mg n systems selected from the SDSS that pro- 
duce strong redden ing signatures in the o ptical spectrum of the 
background quasar. Isirianand et al.l (120081) have discussed two of 
these systems (z a b s = 1.3265 system towards J0850+5159 and 
z abs = 1.3095 towards J0852+3435) that show both the 2175 A 
UV bump and 21-cm absorption. The other two systems that 
show strong signatures of the reddening due to dust but inter- 
estingly are not detected in the 21-cm absorption are discussed 
below. 

1. z abs = 1.3237 system towards J0956+4046: This system is 
amongst o ne of the 12 high-co nfidence 2 175 A absorber can- 
didates of iJiang et all d201 ll) . In the top panel of Fig [T2] 
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we show extinction curve fits to the S PSS spectrum of 
this obj ect using the method described in iNoterdaeme et al.l 
d2009ah . The observed spectrum is well fitted using LMC2 
extinction cur ve and E(B-V) = 0.27 and the bump strength 
(as defined bv lJiang et alj|20l"it) of Ab um p = 0.86. The main 
difference in our procedure is that we have fixed the 
wavelength and wid t h of t he 2 175 A feature as given by 
iFitzpatrick & Massal (120071) . Using all the 491 QSOs with 
Zem within 0.004 to J0956+4046 we confirm the presence 
of bump at the high significance level (see left panels in 
Fig [TZl i. If extinction per hydrogen atom in this system is 
similar to LMC then we expect N(H i) ~ 5 x 10 21 ctrT 2 . 
The estimated E(B-V) value and rest equivalent widths of 
various species measured in this system are identical to 
the 21 -cm absorber at z abs = 1.3265 towards J0850+5159 
(ISrianand et al.|[2008h where the measured integrated 21 -cm 
optical depth is 15.3 km s _1 (see Table 3 of G09). As shown 
in Fig. [TJ J0956+4046 is resolved into two components in 
our GMRT image. The brightest component having a flux 
density of 53 mJy coincides with the optical position. A 21- 
cm absorber like the one we have seen towards J0850+5159 
would have been easily detectable towards it. From Table lB.ll 
we see that the 3cr limit on J rdv is 0.35 km s _1 for this 
component. This together with the above estimated N(H i) 
value suggest that T s /f c > 7800 K. In the case of high-z 
CO absorbers E(B-V)/W(H i) is found to be much higher 
(upto an order of magnitude) c ompared to what is seen in 
LMC dNoterdaeme et al.l2010al Ledoux et al. in prep.). If we 
use the mean E(B-V)/W(H i) found for the CO-DLAs whose 
SED is also well fitted with LMC2 extinction curve we get 
N(H i)~4xl0 20 cirT 2 and T s /f c > 600 K. Therefore, the lack 
of 21 -cm absorption could either be due to the covering fac- 
tor of the gas being much less than 1 or the extinction per 
hydrogen atom in this absorber being much higher than that 
has been seen in LMC. As dusty regions are usually cold and 
the fact that our spectra is sensitive enough to detect the cold 
gas even when f c is as low as 0. 1 favours the second possi- 
bility. The multifrequency VLBA imaging of this source to 
constrain f c and the N(H i) measurement will provide impor- 
tant insights into this dusty Mg n system. 
2. z abs = 0.8620 systems towards J1203+0634: B ased on the 
presen ce of strong Ca n absorption and reddening. IWild et al.l 
(120061) have suggested this system as a possible DLA can- 
didate. We find that the SED can be well fitted with the 
SMC like extinction curve without any need for the pres- 
ence of 2175 A bump (See lower panel in Fig. [12}. The 
equivalent width of Mg n and other species are like the 
other three systems mentioned above. The observed E(B- 
V) for SM C like dust-to- g as rati o corresponds to N(H i)> 
10 22 cirT 2 . iKulkarni et al.l d201 lb have reported a possible 
detection (3.6cr leve l) of 10 /j.m Silicate absorption in this 
system. Interestinglv lKanekar et al.l d2009l) have not detected 
21-cm absorption despite having a very good J rdv limit of 
0.07 km s _1 . Lack of 21-cm detection in this system that has 
a high probability of having very large N(H i) is interesting. 
Unlike J0956+4046, J1203+0634 is a strong radio source 
and all the flux density at arcsecond scales is recovered in the 
milliarcsecond scale VCS images at 2.3 GHz. The strongest 
VLBA component contains 66% of the total flux density. 
Rest of the flux is extended within 170pc. Interestingly in the 
8.4 GHz VLBA image the extended component is undetected 
and the unresolved component seen in the 2.3 GHz VLBA 
image resolves in to a core-jet structure over a LS of < 60 pc. 
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Fig. 12. Fitting the spectral energy distribution of J0956+4046 
(top) and J1203+0634 (bottom). The SDSS s pectrum (black) and 
the un reddened SDSS composite (gray) from lVanden Berk et al.l 
d2001h . The best extinction modelling is obtained for E(B- 
V)=0.27 and a LMC2 extinction law for J0956+4046 (resp. E(B- 
V)=0.34 and SMC law for J 1203 +0634) at the redshift of the ab- 
sorber. The shaded a rea represents Ab um D , the difference between 
the best model using lJiang et al.l (1201 ll) parametrisation with and 
without bump. The orange regions at y=0 indicate the spectral 
range considered in the fitting process. The right panels show 
the distributions of E(B-V) (using SMC-law: red, LMC2: green 
hashed) and Abump obta ined for a control sample of SDSS QSOs 
(ISchneideretal.ll2010h within Az = +0.004 from J0956+4046 
(resp. J1203+0634). The arrows indicate the positions of the 
considered QSO. 

Even if only 10% of the radio emission passes through the 
absorbing cloud the J rdv limit achieved is sensitive enough 
to detect a T-100K gas having N(H i)~10 20 cirT 2 . All this 
suggests that the lack of 21-cm absorption in this system may 
also be related to the extinction per H i atom being much 
higher than the Galaxy or SMC. 



6. Discussion 

6.1. Evolution of 21-cm detection rate andn 2 \ 

The nature of Mg n absorbers and their relationship with the 
galaxies has been a matter of debate. Recently [O n] emission, 
a measure of the star formation rate, has been detected in the 
individual as we ll as stacked SPSS spec t ra of the W T >0.7A 
Mg n absorbers (INoterdaeme et al.ll2010bt iMenardet aDl201ll) . 
The aver age emission li ne lum inosities are found to correlate 
with W T . iMenard et al.l d201 lb show that Mg n absorbers not 
only trace a substantial fraction of the global [O n] luminosity 
density but also recover the overall star formation history of the 
Universe up to z ~ 2. They argue that this is suggestive of the 
scenario where most of the strong Mg n absorbers arise from the 
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Fig. 13. Number of 21 -cm absorbers per unit redshift range, 
«2i(z), for T =0.3 km s . Open circles are for no correction 
for partial coverage. Filled circles correspond to «2i after the 
optical depths have been corrected for partial coverage assum- 
ing f c = Cf. The «2i based on DLAs at z > 2 are plotted as 
squares and crosses with the latter corresponding to values af- 
ter correcting for partial coverage. The curve for non-evolving 
population of 21 -cm absorbers normalized at «2i(z = 1.3) is 
also plotted. Lines and dashed ar eas show the nu mber of ab- 
sorbers per unit redshift for DLAs ([Rao et al. 2006, ; dotted line) 
and strong Mg 11 absorbers (iProchter et al. 2006, ; dashed line). 
Triangles are the number per unit redshift o f DLAs measured 
from the sample of Not erdaeme et al.1 (l2009bl) . 



starbu rst driven galactic outflows (however see, Lopez & Chenl 
120121 who suggest that the W t vs. [O 11] luminosity correlation 
can also be explained by the W r vs. galaxy impact parameter 
anti-correlation and does not necessarily support the 'outflow' 
scenario ). 

Since the volume filling factor of different phases of ISM 
depends sensitively on the stellar feedback, in the 'outflow' sce- 
nario we expect CNM filling factor and the 21 -cm detection rate 
of stron g Mg 11 absorbers to evolve strongl y with the starburst in- 
tensity ( de Avillez & Breitschwerdtl2004l) . But in our analysis of 
the 21 -cm absorption sample presented in the previous sections, 
we find that the 21 -cm detection rate, C, is constant over the 
redshift range of 0.5 < z < 1.5 i.e. 30% of the age of Universe. 

As described in G09, we can estimate the number density per 
unit redshift, n2i(7~o,W ,z), of 21-cm systems with integrated 
optical depth T21 > To and Mg 11 equivalent width W T > W a 
from the number density per unit redshift of Mg 11 absorbers, 
»Mgn(Wo,z), using the equation, 



The H21 for different values of To are giverfl in the last column 
of Table [5] In Fig. Q~3] we also plot the «2i for To =0.3 km s . 
The constancy of C as a function of redshift implies that the 
«2t increases with redshift with a same slope as tiuga- Due to 
the small number of positive 21-cm detections in Mg 11 systems 
the errors on «2i are large and on its own the estimates would 
also be consistent with the non-evolving population of 21-cm 
absorbers (Fig. H3l . We conclude that the current data does not 
show any sign of the evolution of CNM filling factor of strong 
Mg n absorbers and «2i simply follows «Mgii over 0.5 < z < 1.5. 
This is puzzling if the strong Mg 11 absorbers are indeed driven 
by the stellar feedback processes that tend to diminish the CNM 
filling factor. 

6.2. Implications for 21-cm detection rates based on samples 
of DLAs at z > 2 

Unlike strong Mg 11 absorbers that are often favoured to be as- 
sociated with the galactic winds, DLAs are believed to probe 
the gas confined to the disks of galaxies. It is also interesting to 
note that, in contrast to the several Mg 11 host galaxy detections, 
the searches to detec t DLA host galaxies have mostly resulted in 



null detections (e.g. iMoller & Warren 19931 Fvnbo etall 



Rahmani et al. 201 0t~ Fvnbo et alT 201 lh Bouche et al 



2010 



2012; 



Noterdaeme et al.l 120121) . Consequently DLAs are believed to 
correspond to low star-forming galaxies or trace the outer re- 
gions of extended H 1 disks far away from the star-forming re- 
gions (iRahmani et al. 2010[). Given that a large fraction (~40%) 
of strong Mg 11 absorbers are likely to be DLAs, it is interesting 
to view the 21-cm detection rates derived for the Mg 11 systems 
in the context of 21-cm absorption searches based on bonafide 
DLAs. 

Since at z>2, Lya can be observed from the ground based 
facilities the systematic searches of 21-cm absorption in DLAs 
have only been possible at z > 2. A summary of all the 28 DLAs 
searched f or 21-cm absorption at 2.0<z a t>s^3.5 is provided in 
Table 6 of dSrianand et al.ll20T2T) . Using the subsample of DLAs 
towards the quasars with the VLBA 'core' identifications from 
this Table, we estimate the 21-cm detection rate, C = 0.22^ ^ 
for To = 0.3km s _1 . In Fig|5]we plot «2i corresponding to this. 
The method used is exactly similar to the one used for estimat- 
ing «2i from Mg 11 systems in Section l6Tl except that the number 
per unit redshift range of DLAs («dla) is used in the place of 
«M g n in Equation|2] The /?dla requ ired for this purpos e has be en 
me asured from the DLA sample of iNoterdaeme et all (l2009bb . 

iKanekar et al.l (120091) used the detection rate of 21-cm ab- 
sorption in Mg 11 systems and the fraction of Mg 11 systems that 
are DLAs to argue that the 21-cm detection rate in DLAs is 
strongly decreasing with redshift. In the W r >1 A Mg 11 subsam- 
ple with the VLBA 'core' identifications presented in this paper, 
the 2 1 -cm detection rates at 0.5 < z < 1 and 1.1 < z < 1.5 for To 
= 0.3kms _1 and/ c = Cf correction are 0.17 + q ?^ and O-IO^qq re - 
spectively (Table©. From the HST sample of lRao et alj (12006), 
the fraction of strong Mg 11 systems that are DLAs in these red- 
shift bins are 27^% and 38^[ \% respectively. If we assume that 
all the 21-cm detections are DLAs then the fraction of DLAs 
showing 21-cm absorption will be 0.63 and 0.26 respectively for 
the low- and high-z bins. Since all the 21-cm detections need not 
be associated with the DLAs, these fractions should only be con- 



«2i(T 2 i > To, W r > W ,z) = C(z) x n Mgll (W r > W ,z). 



(2) 



7 n 2 i for strong Mg n absorbers at 0.5 < z < 1.5 have been es- 
timated by approximating ziMgii as, ;jmbIi(W ,z) = «o x (1 + z) y with 
n = 0.080!°-°^ and 7 = 1.40 + 0.16 for W ■■ 
I2006h . 



1 A (Prochter et al. 
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sidered as upper limits. The 21 -cm absorption detection rate in 
the DLAs at z > 2 is lower by at most a factor 3 compared to the 
upper limits we obtained using the W r >1 A Mg n absorbers at 
0.5 < z < 1.0. However, at 1.1< z <1.5 the detection rate or the 
upper limit from the Mg n systems are consistent with the 21 -cm 
detection rates from the DLAs at z > 2. This would imply that ei- 
ther the cross section of the 21 -cm absorbing gas amongst DLAs 
has increased from z — 3.5 to z — 0.5 or that the significant frac- 
tion of 21 -cm absorbers at 0.5 < z < 1.0 arise from sub-DLAs. 
The N(H i) measurements of Mg n systems with good 21 -cm 
optical depth limits will be useful to distinguish between these 
possibilities (see also lCurran & Webbll2006( lCurranll2012l) . 



6.3. Caveats and future work 

In previous sections, we presented results from 2 1 -cm absorption 
line searches based on the sample of 85 strong Mg n absorbers 
at 0.5 < z < 1.5. When no correction for the effect of partial 
coverage is incorporated we find that the 21 -cm detection rate 
can be underestimated by as much as factor 2. Using milliarc- 
second scale VLBA maps for 53 quasars with 57 intervening 
absorption systems from this sample, we determined the realis- 
tic 21 -cm optical depths and detection rates. This is the largest 
sample of 21 -cm absorbers assembled at any redshift range, and 
we systematically take into account the effect of varying 21 -cm 
optical depth sensitivity and partial coverage towards different 
quasar sight lines. Despite this the small number statistics still 
play a significant role and results are not without caveats. 

First caveat is related to our assumption that f c can be es- 
timated from the VLBA continuum images at 20-cm and 13- 
cm, whereas the redshifted 21 -cm observations are at the wave- 
lengths longer by a factor 2-4 compared to these. The actual val- 
ues of f c are likely to be different at the observing frequency. It 
is not obvious if the effect of this assumption is to simply in- 
crease the scatter or produce systematic wavelength dependent 
offsets as one also needs to take into account the effect of fi- 
nite resolution of the VLBA images when identifying 'core' and 
determining f c . Given the sensitivity of current receivers and the 
difficulties involved with calibration, the low-frequency milliarc- 
second scale imaging and spectroscopy is possible only for a 
very few strong targets. The VLBI imaging at 20-cm or shorter 
wavelengths will be the optimal way to estimate f c for the large 
samples as presented here. Therefore it will be valuable to deter- 
mine the appropriateness of this method by measuring the filling 
factor and extent of the 21 -cm absorbing gas via VLBI spec- 
troscopy of z < 0.2 21 -cm absorbers that is possible with the 
currently available VLBI receivers. 

Second caveat is related to the fact that the 21 -cm absorption 
searches in the different redshift ranges are based on the samples 
of different atomic or ionic species which complicates the overall 
interpretation. For example, the samples at z > 2 are based on the 
DLAs whereas at z < 2 on the Mg n systems which selects both 
the DLAs and sub-DLAs. Measuring A^(H i) in the Mg n systems 
is important to address this issue. More importantly it is well 
known that the optical spectroscopic surveys such as SDSS are 
performed after the color selection of QSO candidates and are 
biased against the dusty sight lines that are relevant for the cold 
atomic and molecular absorption lines th at can be searched at 
radio wavelengths (e.g.. lCarilii et al.fl998b . Therefore, radio ab- 
sorption line surveys based on the samples of optical absorption 
lines provide a very incomplete view of the evolution of the cold 
atomic and molecular gas in the galaxies. The blind searches 
of radio absorption lines are required to overcome this barrier 



and eliminate the biases due to the pre-selection based on differ- 
ent atomic/ionic species. Such large spectroscopic surveys will 
be possible with the upcoming Square Kilometre Array (SKA) 
Pathfinders such as APERTIF, ASKAP, EVLA and MeerKAT. 

We finish with the following question: what is the relation- 
ship of the 21 -cm absorbers with the galaxies? This link is re- 
quired to establish the direct connection with the nearby galax- 
ies that are detectable in H 1 21 -c m emission. Using 21 -cm emis- 
sion line maps of nearby galaxies. IZwaan etaL 1 (12001 have con- 
cluded that sight-lines with log A^(H i)(cnT 2 )>20.3 occur with a 
median impact parameter, p<7.8 kpc. It is also clear from their 
Fig. 14 that when 10kpc< p <15kpc the probability of hav- 
ing high Af(H i) is roughly between 50% to 70%. Galaxies have 
been f ound within 20 kpc to the DLAs and 21 -cm absorbers at 
z<1.0 dRao et al.ll2003ll201 ll) . The redshift distribution of DLAs 
typically require p<30 kpc to b e consistent w ith the observed lu- 
minosity function of galaxies (ISteidellll995l) . While these stud- 
ies constrain typical impact parameter ranges where 21-cm ab- 
sorption is likely to be detected, they do not tell much about the 
nature of 21-cm absorbers itself due to the lack of information 
on spin temperature. Therefore, observations of 21-cm absorp- 
tion towards the bright radio sources at small impact parameters 
to the gas disks/halos of nearby galaxies are required to deter- 
mine the 21-cm absorption cross section of the galaxies and un- 
derstand the nature of 21-cm absorbers (| Ca rilli & van Gorkoml 
1 1 992l iGupta et al.ll20 1 0t iBorthakur et al.ll20 111) . 

7. Summary 

We have presented results from a systematic search of 21-cm ab- 
sorption in a sample of 17 strong Mg n absorbers at 0.5 < z a b s <1.5 
using the GBT, GMRT and WSRT This resulted in 4 new 21-cm 
detections. The VLBA 20-cm maps of 40 quasars with 42 inter- 
vening Mg n systems are also presented. 

Combining 21-cm absorption measurements f or 50 strong 
Mg i i systems presented in this paper and Gupta et akl 
( 2009|) with the m easurements from literature (lLanel l2000l 



iKanekar et al.ll2.0d9b . we assemble a sample of 85 strong Mg n 
absorbers at 0.5 < z < 1 and 1.1 < z < 1.5. Using the milliarc- 
second scale VLBA maps for the 54 quasars with 57 intervening 
absorption systems from this sample, we estimated the actual 
21-cm optical depths and detection rates after correcting for the 
partial coverage. In our detailed analysis of this 21-cm absorp- 
tion sample, we have taken into account the effect of varying 
21-cm optical depth sensitivity and the partial coverage towards 
different quasar sight lines. Our main results are summarized as 
follow: 

1. The 21-cm detection rate is found to be higher towards the 
quasars with flat/inverted spectral index, a<0. More than 
73% of the 21-cm detections are towards the sources with 
linear size, LS<100pc. This may imply that the character- 
istic size of the absorbing gas is less than lOOpc. The low 
detection rate towards the quasars with the steeper a are due 
to poor optical depth sensitivity towards the 'core' i.e. low 
covering factor, f c of the absorbing gas. 

2. We find that the largest velocity widths (AV>100km s~') 
are mainly seen towards the quasars that show extended ra- 
dio structure at arcsecond scales. However, we do not find 
any correlation between the integrated 21-cm optical depth, 
frdv or the width of 21-cm absorption line, AV with the LS 
measured from the VLBA images. This suggests that the ab- 
sorbing gas is patchy with a typical correlation length ~30- 
lOOpc. 
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3. We find that when no correction for the effect of partial cov- 
erage i.e. as estimated from the milliarcsecond maps is in- 
corporated the 21 -cm detection rate can be underestimated 
by as much as a factor of 2. 

4. The 21-cm detection rate is constant over 0.5 < z < 1.5 i.e. 
30% of the age of Universe. The constancy of C as a function 
of redshift implies that «2i increases with the redshift with a 
same slope as «Mgii- Due to the small number of positive 21- 
cm detections in Mg n systems the errors on «2i are large 
and on its own the estimates would also be consistent with 
a non-evolving population of 21-cm absorbers (Fig.Qj). We 
conclude that the current data do not show any sign of the 
evolution of the CNM filling factor of strong Mg n absorbers 
and «2i simply follows «Mgii over 0.5 < z < 1.5. This is 
intriguing given that the strong Mg n absorbers are believed 
to be driven by the stellar feedback processes that tend to 
diminish the CNM filling factor. 

5 . We find that the 2 1 -cm absorption detection rate in the DLAs 
at z > 2 is at most a factor of 3 less than the upper limits we 
obtained using the Mg n absorbers at 0.5 < z < 1.0. This 
would imply that either the cross section of 21-cm absorbing 
gas amongst DLAs has increased from z - 3.5 to z — 0.5 or 
that a significant fraction of 21-cm absorbers at 0.5 < z < 1 .0 
arise from sub-DLAs. 

6. We find that 75% of 21-cm detections have Mg n dou- 
blet ratio <1.1, W(Mg ii2852)/W(Mg ni2796)>0.27 and 
W(Mg iLl2796)/W(Fe ili2600)<1.47. This confirms our pre- 
vious finding that the probability of detecting 21 -cm absorp- 
tion i s higher in the systems with high N(H i) dGupta et al.l 
120091: ISrianand et all |2012|) . and is consistent with the 
3cr level correlatio n between A/(H i) and TJfc noted by 
ICurran et al. l <Eoiob . 

7. We present a detailed discussion on the non-detection of 21- 
cm absorption in the two dusty absorbers that produce sub- 
stantial reddening to the background QSOs. This could be 
either related to the gas covering factor being low or to the 
extinction per hydrogen atom in these systems being much 
higher than what is seen in the SMC and LMC. We prefer 
the second possibility on the basis of larger extinction per 
hydrogen atom seen in some of the high-z DLAs with CO 
detections. 

We discuss caveats of our survey results. The analysis is 
mainly limited by the small number of 21-cm detections. 21-cm 
detecti ons of the order of SDSS DLA sample dNoterdaeme et al.l 
l2009bl) will be possible with the upcoming SKA pathfinders. 
Blind searches of 21-cm absorption lines with these instruments 
will provide a complete view of the evolution of cold gas in 
galaxies and shed light on the nature of Mg n systems and DLAs, 
and their relationship with stellar feedback processes. 
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Table 4. Results from the VLBA data 



R.A. 


Dec. 


Zabs 


rms 


Id 


S 


r 


e 


a 


b/a 


<P 


Sr 


/VLB A / 


LS 


a 


(J2000) 


(J2000) 




(mjy beam" 


') 


(mJy) 


(mas) 


(deg) 


(mas) 




(deg) 


(mJy) 


Cf 

J 


(pc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 










VLBA parameters at 20-cm 
















01 08 96 RzH? 
Ul Uo ZO.OtJZ 


00 ^7 9/1 060 
UU J / zt-.uuu 


1 ^71 
1 .J / 1U 


I 8 
1 .0 


1 


489 





_ 


6 ^6 

U.JO 


Szt 

U. JT 


78 
/ o 


Q^0 
" JU 


Q7/ 


SS/l 

JJ4 


^8 
U.JO 










2 


162 


16.8 


80 


9Q QO 


1 9 
U. 1Z 


78 
- / o 




S^ 

U.J J 














3 


58 


33.1 


88 


91 "3S 
Z 1 .J J 


9Q 

u.zv 


-J J 


















4 


90 


51.4 


87 


1^61 
1 J.U1 


^Q 

U.J7 


8zt 


















5 


106 


65.3 


79 


91 71 
Z 1 . / 1 


6zt 


6zt 
-OH- 










01 SI SI ^6S^ 
Ul J*t JH.JOJJ 


00 07 9^ 9^1 
-UU U / Zj.Zj I 


1 1 80^ 

1 . 1 OUJ 


~K 
U.J 


1 


250 





_ 


6 8^ 
o.o j 


1 6 
U. ID 


^8 
-Jo 


96zt 
ZOt 


QS 
U.VJ 


^S7 

<.J / 


6S 
-U.OJ 


0? 11 S9 9Q0Q 


4-1 zt OS 97 1SS 

Tit UJ Z / .HJJ 


1 116^ 

1 .ttUJ 


9 

u.z 


1 


11 





_ 


zt 9zt 


00 
u.uu 


68 
DO 


1 9S 

1ZJ 


11/ 
U. 1 1/ 


1 1 ^ 
1 1 J 


68 
u.oo 










2 


3 


13.3 


57 


1 QQ 


Szt 


9^ 
-Z J 




OQ* 
V.yjy 






0? ztO 08 1 7S8 
UZ HU Uo. 1 / Jo 


9^ OQ 1 S 797 
-Zj UV 1 J. / Z / 


1 "3617 
1 . JOt / 


9 T, 
Z.J 


1 


4080 





_ 


Zt 7S 
t. / J 


8S 
U.oJ 


3Q 

J7 


69S6 
OZJO 


8S/ 
U.OJ/ 


81 


9S 
-U.ZJ 










2 


1224 


9.9 


22 


zt 8S 


9Q 

U.Z7 


9J. 




6S 
U.OJ 






09 SQ 98 S1 SS 

UZ J7 ZO.JlJJ 


00 1 Q SQ Q8^ 

-UU IV J7.70J 


1 ^70 
1 .J J / u 


^ 
U. J 


1 


231 





_ 


~K S 1 
J. J 1 


^6 
U.JU 


1 
-1U 


9^S 

Z J J 


Q8 

u.vo 


^^0 
<- JU 


ss 

-U.J J 


07 /t9 ~K1 ^87^ 
U / t-Z J 1 Jo / J 


4-^Q 11 ^S 69Q 

TJ7 f4 JJ.OZV 


1 1 18S 
1 . 1M-0J 


T, 
U.J 


1 


74 





— 


S Q8 

J.70 


1 7 
U. 1 / 


97 
-z / 


1 1 zt 

1 It- 


76/ 

u. / o/ 


6^ 
O J 


^s 

U. JJ 










2 


13 


7.6 


155 


6 6S 
u.u J 


8^ 
u.oj 


n 

u 




6S* 
U.OJ 






07 18 OQ 168^ 


4-^0 06 "30 s^n 
tju uu ju.jjj 


1 ztzl70 

1 .TT / U 


i 
u. j 


1 


105 





_ 


98 

zo. J J 


SO 
U.JU 


8S 

o J 


9^1 

Z J 1 


6S/ 

U.OJ/ 


16/1 

4U4 


1Q 










2 


20 


30.1 





1zt ^1 
l^f. J 1 


00 

u.uu 


8S 
-o J 




z!6* 
U.HO 














3 


24 


28.8 


135 


71 Q9 
/ 1 .vz 


9^ 
U.Z J 


7Q 
- ly 










08 00 ^6 096Q 
Uo UU JO.UZOV 


4-SO 1 11 9Q0 
TJU 1U t-'+.ZVU 


1 ZL 1 zt6 
1 .HI HO 


^ 

U. J 


1 


81 





_ 


zt 99 
4.ZZ 


Ozt 

u.ut 


J J 


1 1 8 
110 


81/ 
U.O 1/ 


986 
ZOO 


08 

u.uo 










2 


14 


33.6 


24 


1 6 1 
ID. 1 


1 Q 
U. 1 y 


1 Q 

1 J 




6Q* 

u.ov 






08 09 /t8 1^9^ 
UO UZ HO.HJZJ 


_i_9Q 17 ^4 91 1 

TZ7 1 / Jt-.Zl 1 


1 ^618 

1 .JU^O 


T, 
U.J 


1 


7 





— 


~K Qzt 


Zt9 
U.tZ 


79 
- /z 


1 QS 

1 yj 


Ozt* 
U.UH 


<- J J 


88 
U.OO 


08 08 ^Q 6666 
UO UO J7.UODO 


4-ZtQ SO ^6 S9Q 

t^7 JU JU.JZ7 


1 zt071 
1 .hu / 1 


9 1 
Z. 1 


1 


482 





_ 


9 0^ 
Z.UJ 


S1 
U.J 1 


S^ 
-J J 


Q^O 
y JU 


S9 
U. JZ 


^-17 
si / 


9S 
-U.ZJ 


08 1 S 1A 1 691 
UO 1J jt-.IDZ^t 


j.'i'I OS 9Q 01 
tJJ UJ ZV.U1U 


8S1 S 
U.OJ 1 J 


9 
U.Z 


1 


25 





_ 


^ Q6 


1 6 
U. 1 


96 
ZO 


JtJ 


07 
u.u / 


^0 

sjU 


SQ 
U. Jy 


08 1 7 1 S189 

UO 1 / lU.JtOZ 


4-9^ S9 9^ Q70 
tzj jz zj.t* /u 


1 ^060 
1 . JUUU 


^ 

U. J 


1 


105 





— 


T, 81 

J.O 1 


1 9 
U. 1 z 


6S 

OJ 


9 1 S 

Zl J 


84/ 
U.Ot/ 


1 99 
1 zz 


1 S 
U. 1 J 










2 


71 


5.92 


-44 


zt 7^ 

L Y. 1 J 


00 

u.uu 


66 
OO 




ZtQ 
U.t7 














3 


5 


14.5 


-89 


6 67 


00 

u.uu 


S ~K 
J J 










08 IS 06 9S0^ 
UO tJ UO.ZJUJ 


4-zt9 S7 1 8 ^Q^ 

TtZ J / 10.J7J 


1 1 1zt7 

1 . 1 14 / 


"3 
u. j 


1 


164 





_ 


9 90 
z.zu 


ss 

U.J J 


J.6 


91 8 
Z 1 o 


7S 
u. / J 


^18 

Slo 


0^ 

U.UJ 


08 SO 19 91^8 
Uo JW HZ.ZHJO 


4-S1 SQ 1 1 6S1 

TJ 1 J7 1 1 .OJH 


1 ^96S 
1 . JZOJ 


9 
u.z 


1 


65 





_ 


1 80 

1 .OU 


1 
u. 1U 


8Q 

-07 


6^ 
o J 


1 0^ 
1 .UJ 


^-1 S 
slj 


09 

u.uz 


08 S9 11 7^Q1 
Uo JZ / Jy 1 


4-^zt ^S ztO S91 
TJ'r JJ tU.JZI 


1 ^OQS 
1 . JUVJ 


i 
U.J 


1 


47 





_ 


9 ^Q 
Z. jy 


67 
U.u / 


6_t 


60 


7/1/ 

U. / T/ 


60 
OU 


^6 
-U.JO 










2 


4 


7.1 


-124 


8 6Q 


00 
u.uu 


S7 

J / 




68 
u.oo 






OQ ^0 ^S 08 S1 
Vy JU JJ.U0J1 


4-16 zLzL 08 S7Q 
Tt-O f+ UO.J ly 


691 6 
u.ozio 


Zt 

u.h 


1 


171 





_ 


1 1 Zt 
1 l.t 


1 1 
U. 1 1 


9 X 


^^6 
J JO 


U.OJ/ 


09 

J Z 


1Q 










2 


40 


13.5 


-154 


1 ^ 6 
1 J.O 


^Q 
U. Jy 


i J 




S 1 * 
U.J 1 






OQ S^ 97 QS6S 
Vy J J Z / .VJOJ 


4-^9 9S S1 S9S 
TJZ ZJ Jl.JZJ 


1 9^79 
1 .Z J / z 


i 
U.J 


1 


78 







zt 1 ^ 
t+. 1 J 


98 
U.Zo 


1 

- 1 


1 ^9 
1 JZ 


71/ 
U. / 1/ 


1 61 
101 


1 1 
U.14 










2 


16 


19.3 





1 1 
1 1 .u 


^zt 
U. J^+ 


1 

- 1 




SQ 

U.J7 






1 07 18 07^ 1 
1U U / lo.U / Jl 


4-99 S1 96 QSS 
tZZ Jl ZO.VJJ 


S609 
U. JOUZ 


9 

u.z 


1 


150 







8 1 6 
o. ID 


S9 
U. JZ 


6/1 


^SO 
J JU 


ss/ 

U.J J/ 


10 


si 

U.J4 










2 


42 


7.6 


-112 


1 7 8 
1 / .o 


ztS 


7 
- / 




o /n 

U.HJ 






1 O S8 1 1 O/l 97 
1U JO 1 J.U4Z / 


4-/1 Q ^Q ^6 O/l 6 
+Vy Jy J0.U40 


1 9 1 90 

l.zlzu 


n 9 
U.z 


1 


5 







1 S 1 
1 J. J 


o oo 

U.UU 


6 1 

04 


9/1/1 
Z44 


O 'ill 
U. J J/ 


8/10 
04U 


O 70 
U. /U 










2 


27 


20.9 


92 


7.98 


0.00 


90 




0.11* 














3 


23 


35.8 


97 


5.26 


0.00 


-81 


















4 


6 


56.9 


-99 


22.7 


0.00 


64 


















5 


20 


66.7 


-107 


12.7 


0.98 


-62 










11 00 21.0333 


+ 16 29 14.661 


0.8540 


0.4 


1 


201 







3.20 


0.41 


9 


265 


0.76 


<25 


0.22 



Table 4. continued. 
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Notes. Columns 1 and 2: right ascension and declination of the 'component- 1' (see column 5) fitted to the quasar, respectively; Column 3: absorption redshift; Column 4: rms in the map in 
mJy beam 1 ; Column 5: component id; Column 6: flux of the component in mjy; Columns 7 and 8: radius and position angle of the component with respect to 'component- 1', respectively; Column 
9: major axis in mas of the Gaussian component fitted to characterise the milliarcsecond scale structure; Columns 10 and 1 1: axial ratio and position angle of the Gaussian component, respectively; 
Column 12: flux density in mJy from FIRST/NVSS (S7-); Column 13: ratio of total flux density detected in VLBA map and S T (called /vlba)- In case of quasars with the multiple Gaussian 
components the c/ i.e. the ratio of the flux density of the strongest component detected in the VLBA map and the St is also given.; Column 14: largest linear size (LS) i.e. separation between the 
farthest components of the quasar in pc at absorber redshift. For the quasars represented by a single Gaussian component we take major axis of the deconvolved component as the upper limit on 
quasar size; Column 15: spectral index, Soc v~ a , between the 20-cm and 50-cm/90-cm using flux densities from the arcsecond scale maps. 
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Appendix A: Log of the 21-cm absorption 
observations. 

Table A.l. Log for the GBT, GMRT and WSRT observations. 



Quasar 


Telescope 


Date 


Time 
Chrs~l 


J0457-2324 


GBT 


2010 Mar 12 


0.8 


J0800+5010 


GMRT 


2009 Jun 05 


4.4 


J08 17+2352 


GMRT 


2009 Jun 23 


2.4 


J0930+4644 


WSRT 


2009 Aug 22 


22 






2009 Aug 23 








2009 Sep 06 








2009 Sep 08 




J0956+4046 


GMRT 


2009 Aug 09 


12 






2009 Aug 24 




J1007+2251 


GBT 


2009 Oct 15 


4.4 




WSRT 


2009 Oct 24 


12 


Jl 148+5254 


GBT 


2009 Oct 16 


4.3 






2009 Nov 19 








2010 Jan 08 




J1216+5843 


GBT 


2009 Nov 19 


1.9 






2010 Feb 18 






WSRT 


2009 Aug 20 


8.8 


J 1252+4427 


GBT 


2010 Jan 19 


2.3 


J1329+1053 


GBT 


2010 Mar 12 


4.1 






2010 Mar 13 








2010 Mar 15 








20 10 Mar 22 






WSRT 


2009 Sep 16 


9.8 


J1333+1649 


GBT 


2010 Jan 09 


2.4 






om n A/Tor i o 
zuiu iviar iz 




J1408-0755 


GMRT 


2009 Jun 21 


2.5 






2010 Jan 22 + 


4.5 


J1410+6141 


GBT 


2009 Sep 16 


5.0 






2009 Sep 17 








2009 Oct 16 








2009 Oct 19 




J1501+5619 


GMRT 


2009 Aug 24 


3.1 


J1636+2112 


GBT 


2009 Sep 25 


2.4 


J2031 + 1219 


GMRT 


2009 Jun 05 


2.8 


J2129-1538 


WSRT 


2006 Jul 12 


27 






2006 Jul 13 








2006 Jul 17 





Notes. The last column gives the total on-source time. Higher reso- 
lution spectrum of the absorber towards J1408-0755 to study the vari- 
ability. 



Appendix B: Determining flux densities and 21-cm 
optical depths 

For the quasars observed with the GB10 the rms confusion due 
to other radio sources in the beam can be the dominant effect 
that limits the accuracy with which flux density of the back- 
ground quasar and therefore the 21-cm optical depth can be de- 
termined. We use interferometric images from our WSRT obser- 
vations and literature to address this issue. First consider the two 
quasars with 21-cm absorption detections from the GBT data. 



Table B.l. Summary of our 21-cm absorption measurements. 
Listed from left to right are the quasar name, flux density in mJy 
at the redshifted 21-cm frequency, spectral resolution in km s _1 , 
spectral rms in mJy beam -1 channeL 1 at the resolution given 
in the previous column, 3<x optical depth limit for the spectra 
smoothed to 10 km s , and the integrated 21-cm optical depth 
or in case of non-detections the 3<x upper limit for a spectral 
resolution 10 km s _1 . 



Quasar 


Flux 


Sv 


rms 


T 3o-,10 


Jrdv 


J0457-2324 


2400 


3.9 


3.2 


0.002 


0.20±0.02 


JUOUUTJU1U 


126 


4 


2 1 


0.033 






244 


J . o 


2 4 


n ni q 


iq 

^o. i y 


J0930+4644 


415 


3.4 


5.3 


0.024 


<0.24 


J0956+4046(P1) 


53 


3.8 


0.9 


0.035 


<0.35 


J0956+4046(P2) 


23 


3.8 


0.9 


0.080 


<0.80 


J1007+2251 


397 


3.2 


1.7 


0.008 


<0.08 


Jl 148+5254 


122 


3.8 


1.3 


0.021 


<0.21 


J1216+5843 


537 


3.6 


2.7 


0.010 


<0.10 


J1252+4427 


497 


3.9 


2.0 


0.008 


2.06±0.13 


J1329+1053 


138 


3.5 


1.2 


0.017 


<0.17 


J1333+1649 


486 


3.6 


2.1 


0.008 


<0.08 


J1408-0755 


744 


3.8 


1.9 


0.005 


0.51±0.05 




685 


1.9 


1.9 


0.004 


0.48+0.04 


J1410+6141 


93 


3.6 


1.5 


0.031 


<0.31 


J1501+5619 


223 


3.8 


2.1 


0.018 


<0.18 


J1636+2112 


456 


3.7 


2.1 


0.009 


<0.09 


J2031+1219 


595 


3.5 


2.7 


0.009 


2.11+0.15 


J2129-1538 


341 


3.4 


3.2 


0.017 


<0.17 



For J1252+4427, in the FIRST catalogue, there are no strong 
(>10mJy) sources within 15' and the total flux density due to 
other sources within the beam FWHM is only ~20mJy. This 
suggests that the confusion is not a serious issue for this quasar, 
and indeed the flux density estimated from our GBT observa- 
tions matches within ~5% with the estimate from the interpola- 
tion of flux densities from the FIRST (at 20-cm) and WENSS 
(at 92-cm) surveys. The other quasar with 21-cm detection i.e. 
J0457-2324 is a blazar and known to exhibit variability (see 
Section l2~2l . We measure the flux density of 2.72 Jy from our 
GBT observations. The contribution due to other sources within 
the beam FWHM at 20-cm based on the NVSS is 384 mJy 
(no beam-correction). The 90% of contribution comes from a 
source that has a steep spectral index of 0.9. Assuming that this 
other source is not variable would imply that the flux density of 
J0457-2324 is 2.4 Jy (beam correction applied) at the redshifted 
21-cm absorption frequency. 

Now consider systems towards the quasars with no 21- 
cm detections from the GBT. The flux density estimates for 
J1007+2251, J1216+5843 and J1329+1053 at the redshifted 21- 
cm frequency are estimated from our WSRT observations which 
took place within 6 months of the GBT observations (Table lATl i. 
For Jl 148+5254 and J1410+6141, the flux densities at the red- 
shifted 21-cm frequency have been estimated by interpolating 
between the flux densities from the FIRST and WENSS sur- 
veys at 20-cm and 92-cm respectively. The flux densities for 
J1333+1649 and J1636+2112 have been determined using the 
FIRST and TXS surveys. Table lETTI lists the flux densities and 
optical depth values for all the Mg n systems observed. 



8 GBT beam FWHM at 800 MHz, corresponding to the central fre- 
quency of PF1-800 receiver, is 15'. The rms confusion corresponding 
to this is 90 mJy. 



Appendix C: Strong Mg n systems selected f rom kanekar et al.l (]2009t> 



Table C.l. Details of the strong Mg n systems {W x >1 A) and 21-cm absorption measurements taken from lKanekar et aO (120091) . Listed from left to right are the quasar names, 
their emission redshifts, redshifts of the intervening Mg n systems, rest frame equivalent widths of the Mg n/12796, Mg n/12803, Mg 1/12852 and Fe H/12600 absorption lines, 
integrated 21-cm optical depth or in case of non-detection the 3<x upper limit to it for a spectral resolution of 10km s _1 , and quasar morphology at 1.4GHz. The quasars are 
classified as compact (C) when deconvolved size in the FIRST is less than 2" and resolved (R) when the deconvolved size > 2". When the FIRST is not available, the classification 
is either based on the radio morphology from the references provided below or on the NVSS/SUMSS when nothing else is available. The last three columns give /vlba, Cf and 
the largest linear size of the quasars estimated from the milliarcsecond scale resolution maps. 



Quasar 


2em 


Zabs 


W r (Mg nA2796) W t (Mg nA2803) 


W r (Mg 02852) 


W r (Fe 1^2600) 


frdv 


Mor. 


fvLBA 




LS 








(A) 


(A) 


(A) 


(A) 


(km s" 1 ) 








(pc) 










Systems at 0.5 < z < 


I 














R/VYSQ A(Y7 


1 A1Q 


U.o4oj 


2.47±0.01 v 


2.34+0.01 


U.JD+U.Ul 


1 on_i_n n i 


<U.4V 


r 1 
i_ 








DAi nO_i_ 1 7A 


1 1 

L. 1 J I 




1.75 s 


1.20 


<u.zu 


1 no 




K 








RO?zLO— 060 

DUZtu UUU 


1 ROS 


1 

U.JO 1U 


1.90±0.02 v 


1.38+0.02 


U.J71U.UJ 


1 70+0 0? 


<rf) 07S 


r 1 


i nn 

1 ,uu 


1 00 


<-7Q 




1.805 


0.7550 


1.67±0.01 v 


1.48+0.01 


0.78+0.01 


1.25±0.04 £1 


<0.061 


C 


1.00 


1.00 


<88 


B0244-128 


2.201 


0.8282 


1.66±0.01 v ' 


1.56+0.01 


0.49+0.01 


1.33+0.01 


<0.13 


R 3 


1.00 


0.92 


162 


B0409-045 


1.684 


0.8797 


1.58+0.22 


1.16+0.21 


<0.23 


1.21+0.20 


<0.14 


C 1 








B0445+097 


2.108 


0.8392 


3.17 s 


2.19 


0.91 


1.97 


<0.054 


R 2 


0.22 


0.17 


80 


B0812+332 


2.426 


0.8518 


2.67+0.28 


2.19+0.31 


0.55+0.15 


1.12+0.27 


<0.11 


c 


0.07 


0.07* 


<30 


B0957+003 


0.905 


0.6722 


1.96+0.16 


1.57+0.17 


0.32+0.16 


1.15+0.24 


<0.06 


R 








B 1012+022 


1.375 


0.7632 


1.53+0.05 


1.28+0.05 


0.33+0.05 


0.71+0.05 


<0.06 


R 








B 1200+068 


2.182 


0.862 


5.29+0.27 


4.88+0.26 


3.04+0.27 


4.08+0.31 


<0.08 


C 


1.00 


0.66* 


170 


B1210+134 


1.139 


0.7717 


1.22+0.09 


1.14+0.09 


0.38+0.10 


0.87+0.10 


<0.052 


C? 


0.20 


0.20 


<74 


B 1222+438 


1.075 


0.7033 


1.01+0.25 


0.70+0.25 


<0.16 


1.22+0.26 


<0.42 


R 


0.70 


0.70 


<72 


B 1324-047 


1.882 


0.7850 


2.62+0.02 1 ' 


2.36+0.01 


0.78+0.02 


1.87+0.02 


<0.48 


C? 


0.38 


0.38 


<75 


B1343+386 


1.852 


0.8076 


1.61+0.12 


1.45+0.12 


<0.14 


1.20+0.12 


<0.06 


R 


0.70 


0.30* 


81 


B 1402-012 


2.518 


0.8901 


1.14+0.09 


1.03+0.09 


0.30+0.10 


1.08+0.09 


<0.09 


C 


0.65 


0.65 


<78 


B1611+343 


1.397 


0.6672 


1.21+0.06 


1.09+0.06 


<0.07 


0.58+0.07 


<0.09 


C 


1.00 


1.00 


<70 


B1701+593 


1.798 


0.7238 


1.67+0.22 


1.17+0.22 


0.80+0.20 


1.24+0.24 


<0.10 


R 
















Systems at 1,1 < z < 1,5 














B 1005-333 


1.837 


1.3734 


0.93 £2 






0.84 


<0.11 


c 


1.00 


0.78 


280 


B 1136+408 


2.366 


1.3702 


1.37+0.21 


0.69+0.23 


<0.23 


0.60+0.23 


<0.13 


c 


1.00 


1.00? 


102 


B1142+052 t 


1.345 


1.3431 


2.15+0.11 


1.61+0.12 


1.29+0.14 


1.40+0.17 


0.557 


R 


0.60 


0.32 


254 


B 1204+399 


1.518 


1.3254 


1.41+0.11 


1.10+0.11 


<0.11 


0.76+0.14 


<0.23 


C 


1.00 


0.72* 


564 


B 2003 -025 


1.457 


1.2116 


2.65+0. 14 A 


2.17+0.15 


<0.31 


1.27+0.14 


<0.022 


R 









Notes . (T) This system is commo n with G09. (v) Equiv alent widths measured using VLT archival spectra. (S) Equivalent widths from Ba rthel et al.l (Il990f). (r) Equivalent widths from lYork et al.l 
d2007h . <A) Equivalent w i dths f rom lAldcroft et ail i 19941) . <£l) Fe H/12600 is not covered in the VLT UVES spectr um, the HUFe n/t26 00) is taken from lEllison et al.l j2004h . (£2) Equivalent w idths 
taken from lEllison et al.l j2004h . (1) Compact on the basis of NVSS/SUMSS. (2) Projected linear size >2" at 5 GHz jBarthel et aOl988h . <3) Projected linear size >2" at 5 GHz jReid et al.ll 19991) . 



Appendix D: Strong Mg n systems selected from Lanel (|2000h 



The Table 6 of G09 lists strong Mg n systems with 21-cm me asurements from Lane (2000). Only systems at z a b s >0.5 are considered here (see Section |2~3V The systems towards 
B0109+176 and B0957+003 are common with the sample of iKanekar et al.1 (1 2009) and are listed in the Table IC.ll with the more sensit ive optical depth limits. In addition to the 
absorber towards B 1622+238, we have excluded the system towards the quasar B0827+243 as it is non-Mg n selected (see lLane1l2000l for the details). 



Table D.l. Sample of Mg n systems with W r >1A from lLanel d2000) . The column listings are same as in Table IC.l" 



Source 


-■urn 


Zabs 


W r (Mg H/12796) 


W r (Mg H/12803) 


W,(Mg 1/12852) 


TT T /T -1 f\f\\ 

VK r (Fe H/12600) 


jTdv 


Mor. 


fvLBA 


c f 


LS 








(A) 


(A) 


(A) 


(A) 


(kms" 1 ) 








(PC) 


B0109+200 


0.746 


0.5346 


2.26 


1.71 






<0.27 


C 








B0229+341 


1.240 


0.7754 


1.92 


2.02 


<1.13 




<0.43 


R 3 








B0235+164 


0.940 


0.5238 


2.42 


2.34 


0.91 


1.79 


13.0 f 


R 4 


0.3-1 


0.3-1 


<77 


B0420-014 


0.915 


0.6330 


1.02 


0.86 


<0.36 




<1.59 


C 5 








B0454+039 


1.343 


0.8597 


1.53 


1.40 


0.37 


1.11 


<0.11 


C 


0.88 


0.55 


122 


B0805+046 


2.876 


0.9598 


1.01 


0.83 


<0.60 


0.24 


<0.71 


R : 








B1218+339 


1.519 


0.7423 


1.34 


1.08 


<0.70 


1.00 


<0.20 


R 


0.07 


0.07 


<73 


B 1327-206 


1.169 


0.8530 


2.11 


1.48 


<0.40 


0.76 


<2.04 


c 








B1354+258 


2.004 


0.8585 


1.00 


0.86 


<0.10 


<0.20 


<1.05 


R 








B 1556-245 


2.815 


0.7713 


2.07 


1.91 


1.07 


<0.20 


<0.68 


C" 








B 1629+ 120 


1.795 


0.5313 


1.40 


1.35 


0.31 


0.70 


0.49* 


c 


0.30 


0.10 


7165 7 






0.9004 


1.06 


0.67 


0.44 


0.63 


<0.82 


c 








B2212-299 


2.703 


0.6329 


1.26 


1.00 


0.36 




<2.07 


c 1 









(tt 



Notes. (T) 21-cm absorption data are from IWolfe et al.l(ll978h 
at 5 GHz jBarthel et al.ll 1988b . (3) 

(6) Compact in the 5 GHz images from lReid etal] jl999h and lBarthel eTall hoOOh . (7) 



21-cm absorption data are from IKanekar & Chengaiuil (120031) 
.1 119881) . (j) Based on lSwarup et all dl984l). (4) Based on the 1.64GHz image fromlMurphv et all Jl993h . 



in 

(5) 



dUlvestad etal .1119811) 
subarcsecond scale image at 408 MHz. 



Compact on the basis of NVSS. ' Projected linear size >2" 
Unresolved in the 3" scale resolution VLA image at 1.4 GHz 



Based on 2.3 GHz image from[ballacasa et al. ( 1998); see Sai kia et al.1 d 19901) for the 
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Appendix E: 20 cm VLBA maps of quasars from our sample 
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Fig.E.l. Contour plots of the VLBA images at 20cm. The rms in the images are listed in Tableland the maps are centred at the 
Gaussian 'component- 1' fitted to characterise the milliarc second scale structure. The restoring beam, shown as the ellipse, and the 
first contour level (CL) in mJy beam -1 are provided at the bottom of each image. The contour levels are plotted as CLx(-l, 1, 2, 4, 
8,...)mJybeam . 
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Fig.E.l. Continued. 



